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In Abrasive Jet Machining (AJM) fine-grained abrasive 
particles mixed in suitable proportion with high pressure 
carrier gas f usually air, are directed through a suitably 
designed nozzle. The nozzle imparts high velocity to the 
carrier gas and thereby to the abrasive particles at the 
nozzle exit. This high velocity air-abrasive mixture jet 
is then directed onto the workpiece surface. Material remo- 
val occurs due to erosion caused by the impact of high velo- 
city abrasive particles onto the workpiece. 
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Literature survey reveals that erosion caused by 
impact of solid particles entrained in a fluid stream is 
a complex phenomenon and is not fully understood. Further, 
adequate data of its governing parameters are not available. 
Some of these parameters are interdependent and difficult 
to control. The erosion study in such situations in volves 
determination, from the fluid flow conditions, of the number, 
direction and velocity of abrasive particles striking onto 
the workpiece surface and then the estimation of material 
removed from the workpiece. In the present work an attempt 
has been made to carryout theoretical as well as exp er imental 
study of the erosion phenomena in UM. 

In the first phase of the present study, the flow 
governing equations in dimensionless form for two-phase 
flow in the nozzle as well as jet region have been developed. 
These non-linear differential and algebraic equations have 
been solved to estimate the velocity of abrasive particles 
in the air-abrasive jet.J^ 

During solid particle impacts , erosion loss in brittle 
materials occurs due to intersection of lateral cracks from 
adjacent impacts. These cracks are produced on the work- 
piece below the impact site due to the residual stress field. 
Considering a quasi-static approach and using semi -empirical 
fracture mechanics relationships, the expressions for lateral 
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crack length (C^) and. lateral crack depth (z-^) have "been 
obtained in terms of various impact parameters in the second- 
phase of the work. Subsequently, a theoretical model of 
material removal rate has been developed. 

In polishing, deourring and other finishing operations 
the volume of eroded material is of importance , while in micro 
drilling and cutting operations it is the erosion depth 
(depth of penetration) which is more relevant. In the third 
phase of the study, experimental investigations have been 
carried out to investigate the effects of input parameters 
such as grain size, mixture ratio, nozzle pressure, stand- 
off distance etc. on tho volumetric material removal, rate, 
penetration rate, eroded cavity diameter etc. For this 
purpose an experimental Abrasive Jet Machine has been 
designed and fabricated which provides for precise control 
of various input parameters. 

The performance characteristics of machining operations 
in ATM also depend on the nozzle geometry. Tho nozzle impart 
high velocity to the expanding carrier fluid (air) and there 
by to the abrasive particles. High velocity air- -abrasive 
particle mixture while expanding causes wear at the inner 
stir f ac e of the nozzle. Nozzle wear is , therefore , of interest 
not only from the point of view of cost but also from the 
point of view of characteristic performance of the unit and 
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machining accuracy. In the fourth phase of the study, 
experiments have been carried out to investigate the 
effects of nozzle diameter, nozzle length and entrance 
angle on nozzle wear. Other parameters which have been 
considered are mixture ratio and abrasive particle size. 
Experiments have also been carried out with split nozzles 
to examine the wear profile of the inside nozzle surface. 

Experimental results indicate that both volumetric 
material removal rate and penetration rate depend on stand- 
off distance, mixture ratio, nozzle pressure and grain size. 
Both volumetric material removal rate as well as penetration 
rate first increase with increase in stand-off distance and 
then decrease giving an optimum value. Theoretical curves 
also indicate this and there appears qualitative agreement 
between the two. The experimental results further reveal 
that the maximum values of material removal rate and pene- 
tration rate occur at different values of stand-off distance 
with maximum material removal rate occurring at a higher 
value. This indicates that operations such as deburring 
and finishing should be carried -out at larger stand-off 
distances, whereas smaller stand-off distances should be 
chosen for micro -drilling etc. 

Experiments on nozzle wear indicate that wear increa- 
ses when nozzle length is increased and when nozzle diameter 
and entrance angle are decreased. Increase in mixture ratio 
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and size of the abrasive particles also increase the 
wear rate. 

In the concluding portion of the work, the theoretical 
as well as experimental results are discussed in detail in 
order to understand the basic mechanics involved in the 
process . 



CHAPTER 1 


INTRODUCTION AND REVIEW OP IITERATDRE 
1 . 1 Intro duct ion 

Recent advances in science and technologsr specially 
in the areas of nuclear and aerospace engineering have imposed 
a severe demand on the material requirements. As a result 
of this, the metallurgists and material scientists have been 
continuously developing newer and larger varieties of mate- 
rials. These materials are of increasingly higher strength 
and are harder and difficult to machine. The advent of such 
high strength, temperature resistant, hard and brittle mate- 
rials posed the serious problems to the manufacturing engineers. 
Besides, the processing of the parts with higher accuracies 
and complicated shapes have also been difficult, time consuming 
and uneconomical bsr conventional techniques of machining. 

In view of the seriousness of these problems the need for 
developing newer concepts in machining was emphasized [l]. 

To meet these challenges several new manufacturing 
processes have been developed during the last two decade 
or so. Unlike the conventional machining processes these 
processes are not affected bjr hardness, toughness or brittle- 
ness of the material. Purthcr, these processes can produce 
intricate shapes with high degree of accuracy on any workpiece 
material through suitable control of the various physical 
parameters of the process. The newer machining processes 
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so developed are often called 'Unconventional Machining 
Processes' or 'Modern Machining Processes'. These are 
unconventional in the sense that they do not employ conven- 
tional tools for material removal , instead energy in its 
direct form is utilized. 

The underlying principle of these modern machining 
processes is to apoly some sort of energy, namely mechanical, 
electrical, chemical, thermal, or magnetic, to the workpiece 
directly and have the desired shape transformation or mate- 
rial removal from the work surface through use of known 
scientific principles. One of these new technological 
processes involves material removal by erosion or abrasion 
due to surface impact of hard particles (abrasives). The 
erosion of work material by impact of abrasive particles 
necessitates pneumatic or hydraulic pressure and fluid 
stream for carrying these abrasive particles. The process 
employing this mechanism of material removal is the Abrasive 
J et Mach in ing ( AJM ) . 

1.1.1 Abr as ive J et Machining 

Abrasive jet machining is a process in which material 
removal takes place due to the erosive action of a stream of 
fine grained abrasive particles impacting at high velocity 
onto the work surface. Pine grained abrasive particles are 
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fed in suitable proportion from an abrasive feeder into 
the mixing chamber where it is mixed with high pressure 
gas or air (from compressor). The mixture of air and 
abrasive is directed through a pipe to a suitably designed 
nozzle. The nozzle and the nozzle holder are attached 
t,o a nozzle food device to adjust and locate the nozzle 
with reference to the work surface to bo machined. The 
nozzle tip is held at a certain distance from the workpiece 
and in a definite orientation, generally normal to the 
surface to be machined. The distance between the nozzle 
tip and the workpiece surface is generally called the 
stand-off distance. The nozzle used imparts high velocity 
to the expanding carrier fluid and thereby to the jet of 
abrasive particles at the nozzle exit. This high velocity 
air-abrasive jet is directed onto the work surface to be 
machined. Material removal occurs due to erosion caused 
by the impact of high velocity abrasive particles onto the 
work surface. The schematic diagram of abrasive jet machin- 
ing is shown in Tig. 1.1. 

1,1.2 Applications 

Abrasive jet machining process finds its successful 
application for machining of hard and brittle materials. 

The process can also be used for machining of complicated 
shapes and contours. Due to absence of heat during machining, 
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this process has been found particularly suitable for 
machining of heat sensitive materials, lavoie [2] and 
Ingulli [3] have mentioned a number of interesting appli- 
cations of this process. Some of these are summarised 
below. 

(1) Clearing corrosion, oxides, dirt and other contami- 
nants from such diverse type of parts as electronic 
components and leather artifacts. 

(2) Trimming of resistors to specified values. 

(3) Cutting and micro -drilling of thin sectioned fragile 
components made of glass, germanium, refractories, 
ceramics, mica etc. 

(4) Etching of trade marks, identifying numbers or other 
markings on machined items. 

(5) Polishing ox plastic, nylon and teflon components. 

(6) Removing of parting lines and flashings from 

moulded parts and cleaning of intricate mould cavities. 

(7) Debarring holes which are virtually impossible to 
achieve by conventional deburring tools. 

(8) Prosting of interior surfaces of glass-tubes. 

(9) Reproducing of designs on a glass surface with the 
help of masks made of rubber, copper etc. 

(10) Cutting of very hard materials such as diamond etc. 
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1.1.5 Limitations 

The application of the process is limited hy the 
following considerations; 

(1) The capability of the process is limited because of 
low material removal rate (material removal rate for 
machining of glass has been reported to be approximately 
40 gm/min. or 0.0164 cm 3 /min. [4]). 

(2) The tapering effect is inherent characteristic of the 
process because of unavoidable flaring of the abrasive 
jet . 

(3) The accuracy of the process is also affected '.by stray 
cutting which again is difficult to avoid. 

(4) Abrasive may get embedded in the woi'k surface. 

1.2 Variables in Abrasive Jet Machining 

The variables which influence the material removal 
rate, surface finish and accuracy of machining in AJM have 
been summarized below. 

(a) Carrier fluid : 

(i) pressure, 

(ii) viscosity and 

(iii) density, 
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(b) Abrasive : 

(i) type of abrasive grain and 

(ii) size of abrasive grain, 

(c) Nozzle ; 

(i) geometry of the nozzle (i.e. length, diameter, 
entrance angle etc.) and 

(ii) nozzle material, 

(d) Machining variables : 

(i) stand-off distance, 

( ii ) shape of cut , 

(iii) angle of impingement, 

(iv) mixture ratio and 

(v) jet velocity, 

(e) Work material. 

1.2.1 Carrier Fluid 

The carrier fluid used in AJM should be non-toxic, 
cheaply available and capable of being dried and cleaned 
without any difficulty. For better accuracy and efficient 
machining it should be of less flaring characteristic 
after discharging from the nozzle exit. The carrier fluid 
which can be used are air, carbon dioxide, or nitrogen. 
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The air is most widely used because of its cost and easy 
availability. Even clean and dry shop air, directly from 
the air lines, can be used after suitable filtering. 

The pressure of the carrier fluid influences the 
material removal rate in AJM. For machining of glass with 
aluminium oxide abrasive Bhattacharyya[4] has found that 
the material removal rate first increases linearly and 
then tends to saturate with further increase in pressure 
(Fig. 1.2). A similar effect has been reported by Sarkar 
and Pandey [5]. In the existing literature on the subject, 
no mention has been made about the influence of viscosity 
and density of the carrier fluid on material removal rate. 

It is, however, believed that a low viscous and lighter 
fluid would give higher material removal rate in comparison 
with heavy gas or fluid. 

1.2.2 Abras ives 

The choice for abrasive type depends upon the type 
of machining (roughing or finishing) , cost of machining, 
properties of workpiece material, type of operation etc. 

The abrasive grain should be sharp and irregular for bettor 
cutting action and fine enough to remain in suspension with 
the carrier fluid. For cutting purposes aluminium oxide 
and silicon carbide abrasives are generally used. For 
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deburring, polishing and cleaning purposes dolomite, glass 
beads etc. are most suitable. Sodium bi-carbonate being 
hygroscopic is rarely used for light finishing operations. 

The working temperature with sodium bi-carbonate abrasive 
should be kept below 50°C above which it will loose water 
and the exit velocity from the nozzle will be reduced. 

Re-use of abrasive is not recommended because of decreased 
cutting ability and contamination with the abraded particles 
of the work material may cause clogging of nozzle passage. 

Material removal rate in AJM depends upon the 
size and shape of the abrasive grain. Coarse grain being 
more irregular have better cutting action compared to fine^ 
grain abrasive. Also, fine grain abrasives having sticking 
characteristic make the flow irregular and increase the 
choking tendency in the nozzle and flow system. The most 
favourable grain sizes for AJM are 10 to 50 pm. Neema 
and Pandey [6] have shown that for an specified machining 
condition there exists an optimum grain size for maximum 
material removal rate. Their experimental results for 
machining of glass with silicon carbide abrasive are shown 
in Pig. 1.5 which indicates that the optimum grain size 
under these conditions is 40 pm. Pine grain abrasives give 
better surface finish and are recommended for polishing, 
deburring and other finishing operations. Table 1.1 gives 
typical applications of some commercially available abrasives. 



9 


1.2.3 Nozzle 

The abrasive particles mixed with carrier fluid at 
high pressure is allowed to expand through the nozzle. 

During expansion the carrier fluid gains high velocity and 
interacts with the abrasive particles. As a result of 
this abrasive particles also accelerate and attain high 
velocity at the nozzle exit. The fluid particles, having 
smaller inertia compared to abrasive particles, accelerate 
faster. Hence, large length and small entrance angle of 
the nozzle enables the carrier fluid to have sufficient 
time to impart momentum to the abrasive particle. But, on 
the other hand, a larger length will have increased wall 
friction. This will cause further reduction in the particle 
velocity. Thus, there exists an optimum length of nozzle 
which will give maximum particle velocity for an specified 
enirnace angle and loading ratio (the ratio of the mass 
flow rate of particle to the mass flow rate of air). The 
optimum nozzle length for parallel side nozzle of circular 
cross-section with 250 pm silicon carbide particles has 
been found to be of the order of 20 times the nozzle dia- 
meter (Tig. 1.4). Wolak et al. [7], however, found that 
the optimum length to bo independent on the loading ratio 
(fig. 1.5)* 
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The nozzle used in AJM has to withstand severe 
erosive action due to flow of air and abrasive mixture. 
Thus, it must be made of high wear resistant material. 

The common materials for nozzle are sapphire and tungsten 
carbide. 

Depending on the requirements, the nozzle may be 
either straight or right angled with various sizes of 
round and rectangular cross-sections. The dimensions of 
various kinds of round and rectangular slot nozzles of 
different materials and their nozzle life are indicated 
in Table 1.2. 

1.2.4 Machining Variables 

Ingulli [3]> Bhattacharyya[4] and Pandey et al. 

[5 and 6] have shown that the stand-off distance has 
considerable effect on the material removal rate as well 
as the accuracy of machining in AJM. They have found that 
the material removal rate first increases with increase 
in stand-off distance and then decreases with further 
increase giving an optimum value (Fig. 1.6). A large 
stand-off distance results in the flaring up of the jet 
which leads to poor accuracy [3-5]. 
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During the study of erosion of a surface by impact 
of solid particles Wellinger and,- Uet !z . [8]., Finnie [9] and 
Finnie et al. [10 ] have revealed strong dependence of 
erosion rate on the angle of impingement, hardness of the 
target surface and velocity of impacting particles. Sheldon 
and Finnie [ll] in their study have shown that the "brittle 
materials give maximum volume removal rate when the impinge- 
ment angle is nearly 90 degrees, while ductile materials 
give maximum erosion when impingement angle is between 10 
to 30 degrees (Fig. 1.7). 


The kinetic energy of the impinging particles of 
the jet, is utilized for material removal in AJM, which 
is a function of the particle velocity. Thus, material 
removal rate is dependent on the pressure of the carrier 
fluid at nozzle inlet, nozzle design, abrasive type, grain 
size and mixture ratio. ITot only the flow characteristics 
of gas-particle flow system but the number of particles 
taking part into erosive cutting of abrasive jet also 
depends upon the mixture ratio. The mixture ratio a is 
generally defined as the ratio of abrasive particles mass 
flow rate M to the mass flow rate of gas (air) and abra- 
sive particles mixture (M + M ). Thus, 

o Jr 




a 


d.i) 
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At a particular pressure .material removal rate increases 
with increase in abrasive flow rate and is influenced by 
the size of the abrasive particle. After reaching an 
optimum value, the material removal rate decreases with 
further increase in abrasive flow rate [4]. 

1.2.5 Work Mat erial 

Abrasive jet machining is recommended for machining 
of brittle materials such as glass, mica, ceramics and 
refractories. Most of the ductile materials are difficult 
to machine by this process. The material removal rate 
appears to be related to the hardness of the work material 
(Fig. 1.8). 

1.3 Impact Erosion 

It is important to consider the response of target 
material under the action of solid particle impact before 
analysing any erosion process. This response may be 
’ductile* or ’brittle'. In ductile materials, large plastic 
strain precedes fracture and material is removed by the 
displacement action of the impacting particle. By contrast, 

brittle material fractures without observable plastic 

A' 

deformation and material is removed by the propagation 
and intersection of different systems of cracks formed 
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at the impact site. Although the concepts of ideally 
ductile and ideally brittle behaviour are over simpli- 
fication, they do describe to a certain degree of approxi- 
mation the behaviour of many real materials and allow 
analytical solutions to be developed. The mechanisms by 
which ductile and brittle solids erode must be understood 
before one can make realistic analysis for some real 
engineering materials having intermediate type of behaviour. 

1.5*1 Ductile Erosion 

Various experimental and theoretical studies have 
been carried out to understand the process of erosion of 
a ductile metal. Experimental studies on fly~ash erosion 
by Fisher and Davis [12] have revealed that erosion loss 
increases with increase in both impact velocity and stand- 
off distance. Stoker [13] has shown that erosion rate is 
dependent on the impingement angle of the impacting particle. 

Perhaps the first scientific approach for the study 
of erosion phenomenon in ductile materials under the action 
of a stream of solid particles has been made by Finnie [14]. 
The theory assumes that a hard, angular particle, impinging 
upon a smooth surface at an angle 6 from the work surface 
will cut into the surface much like a sharp tool. He 
solved equations of motion of cutting tip of the particle 
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with, the assimiption that the ratio R, of the vertical 
component to the horizontal component of the force acting 
on the cutting tip of the particle, remains constant 
through.- out the impact . Depending on the impingement 
angle ©, two kinds of cutting processes have "been identi- 
fied ; (a) the particle cuts into the ductile surface, 
and subsequently leaves it removing the fragment of chip; 
and (b) the particle is stopped during its scooping action 
at some depth when its kinetic energy is exhausted. The 
limiting value of the impingement angle G Q distinguishing 
between these two situations has been expressed as 

G 0 = tan” 1 ( | ) • (1.2) 


Pinnie [14] has shown that the eroded volume Y from the 

target surface by impacting particles of mass M and 

P 

velocity u . is given by 


Y 



lm 


(1.3) 


where a ^ is the flow stress of the target material and 
f(0) is a function of impingement angle and has been 
expressed as 
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f(e) = Lm- 2 £ - ^Jjti£?§l , 0 <e<e o (1.4) 

and 

f (o) = ( ) . e o < © < 90° (1.5) 

This mode], is quite successful in explaining many 
features of solid particle erosion. However, quantita- 
tive discrepancies appear concerning the effect of flow 
stress, the velocity exponent and the applicability of 
the model itself for impingement angles close to 90 degrees. 
The theory predicts no erosion at normal impacts (© = 90°), 
while experimentally sizable erosion values are obtained. 

In the later work, Finnie and McFadden [15] have extended 
Finnie's original analysis by modifying the interaction 
forces between the particle and the target surface. They 
obtained the velocity exponents of about 2*5 which is close 
to the value found experimentally. 

The erosion rate not only depends on the impingement 
angle but also on the shape of the erosive particles. The 
spherical particles give lower erosion rate than the angular 
particles of the same materia]-. This was verified experi- 
mentally by using different Sand grains [9], 
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The erosion resistance of various metals when 
impacted upon by angular silicon carbide particles of 
250 ).un average size show a linear increase with Viewer's 
hardness on log-log representation (?ig. 1.9) • It has 
also been observed that increase in hardness either by 
cold wording or through heat treatment appears to be in- 
effective in raising erosion resistance [10]. 


By incorporating the idea of particle energy/ - 
threshold, Bitter [16,17] has developed an expression to 
give a better account for erosion at various impingement 
angles. He has postulated two types of wear, ’deformation 
wear* and ’cutting wear', the former is caused by repeated 
deformation, while the later is due to the cutting action 
of free moving particles during collision. Brittle 
materials generally erode due to deformation wear, while 
ductile materials erode due to cutting wear* His approach 
is based on the computation of plastic energy dissipation, 
from the impact parameters of a single erosive particle 
under elastic-plastic collision. He has shown that the 
deformation wear volume Y^ can be expressed as 


M (u . sin ■© 

Y 


u ) 1 

pe y 


2 e 


d 


9 


( 1 . 6 ) 
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where u^ 0 is the elastic impact velocity at which the 
elastic limit is just reached in the target material and 
the deformation wear factor is a material property 
defined as the amount of energy needed to remove a unit 
volume of material. 


The cutting wear exists when particles strike a 
body at an acute angle scratching out some metal from the 
surface. The particle impact velocity u , m is resolved 
into two components, one normal to the target surf an e and 
other parallel to it. The normal component of particle 
impact velocity u^^ causes penetration of particle into 
the target surface, while tangential component of particle 
impact velocit]? - u ^ gives the scratching action. Similar 
to the earlier erosive cutting theory [14], two possibili- 
ties exist for the cutting action. If 0 < © , the incident 
particle leaves the surface with a non-zero velocity after 
cutting, but for © > © 0 the particle velocity is zero as 
its kinetic energy is exhausted, lor the first case 
(© < © Q ), if uh.^ is the residual tangential component 
of the particle impact velocity, the energy absorbed in 
the process of scratching is y M (u 2 . cos 2 © - u'? j.). 

In the second case (© > 9 ), the energy absorbed in the 
process of scratching is 4- cos 2 ©. The cutting wear 

factor s , defined as the quantity of energy needed to 
c 
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scratch.' a unit volume of material from the target surface, 

depends on the mechanical properties of the target material. 

Bitter [17] has shown that the resulting cutting wear volume 

V can he expressed as 
c 
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and 
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( 1 . 8 ) 


Here c^ and c 2 are constants which depend on the physical 
and mechanical properties of the target and the erosive 
particle materials. The eroded volume V due to combination 
of both deformation wear as well as cutting wear can he 
expressed as 


Y 


V. + V 
d c 


(1.9) 
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Hells on and Gilchrist [18] have simplified Bitter's 
model and have given a simpler form of cutting wear equation 
while retaining the deformation wear equation given in 
reference [16]. For different impingement angles they have 
expressed the erosion volume Y as 


Y = 


V d + 


M P (“pin 008 e " “pint) 
2 £ 


; e < ©, 


( 1 . 10 ) 


and 


M p u pim cos 2 © 

v = Y d + * 6 > 6 o > 


( 1 . 11 ) 


where Y^ is the deformation wear volume given by equation 

(1.6) end e is the cutt ing wear factor [17]. 
o 

Tilly [19] has proposed a two-stage mechanism of 
erosion recognizing explicitly that particles impacting 
at near normal incidence may fragment and the fragments 
may subsequently erode the exposed surfaces. He explained 
the reported decrease of erosion with decreasing particle 
size and introduced the concept of particle size threshold 
for effective erosion. The studies on particle size depend- 
ence of erosion [20-22] have further confirmed that erosion 
rate increases with increase in particle size beyond the 
threshold value and ultimately reaches a saturation point. 
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The saturation point, however, is different for different 
materials. Prom the consideration of deformation heat 
generation and its diffusion rate into the target material, 
Shewmon [23] has also shoxm that a particle size threshold 
exists below which no erosion occurs. 

Sheldon and Eanhere [24] have examined the mechanism 
of single particle erosion of ductile materials and have 
developed a method to describe the deformation and machining 
actions observed using indentation theory and energy balance 
equation. Their results give the erosion volume V as 


V = 


c d 3 3 , £P 

c 3 p pirn ^ H* 


3/2 


( 1 . 12 ) 


where c^ is a constant, is the Ticker's hardness value 

of target material, d is the diameter and P is the mate- 

P p 

rial, density of the impacting particle . They obtained 
the velocity exponent as 3 which is notably higher than the 
value of 2 obtained in earlier studies [9, 14-18]. Their 
value of the exponent is, however, comparable to the value 
published for multiple particle erosion [17,22], Earlier 
studies [21] have indicated that the velocity exponent is 
independent of the type of material while later studies [25] 
have shown some dependence on attack angle 6. 
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Hutchings and; Jf inter .[26, 2-7] hay©?, studied the 0 ©ro6ion 
process using large (3 mm diameter) spherical as well as, 
angular particles with particular emphasis on the geometry 
and mechanisms of material removal. The characteristic 
deformation pattern resulting from particle impact consisted 
of a depression and a lip or rim of displaced material. 
Further, they found evidence for a critical particle velo- 
city above which material is displaced from the crater lip. 
The spherical particles cause ploughing while the angular 
particles cause cutting. The particle rake angle at the 
contact was found to be the controlling parameter for cutt- 
ing [28-31]. 

Ives and Ruff [32] have observed that impact at low 
angle of impingement (around 20 degrees) deforms and dis- 
places the material from the crater into a lip at the exit 
and sides, whereas normal impingement produces a more 
uniform lip of material around the crater with considerable 
plastic flow. The material removal under normal impacts 
occurs due to ploughing and flaking [33] followed by fatigue 
[34] under repeated impacts. 

Follansbee et al. [35] have proposed an accumulated 
damage model for the low velocity erosion of ductile metals 
due to multiple impacts of spherical particles at -nor mal 
incidence using an elasto-plastic finite element code. 
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Sundararajan [56] lias also proposed a theoretical model 
for the extent of lip formation as well as the erosion 
loss during single particle impact. He has employed the 
criterion of critical plastic strain at which the defor- 
mation in target material localizes and results in lip 
formation, later Sundararaj an and Shewmoh [5-7] have given a 
new model for erosion of metals. They have obtained the 
velocity exponent as 2,5 for normal impacts. 

levy [58] has observed that the micro structure of 
target material plays an important role in the solid particle 
erosion of ductile alloys. Roley and levy [59] have 
reported that the erosion resistance of steel increases 
with increase in ductility upto a certain limit depending 
on the strain hardening coefficient. They have further 
reported that erosion rate increases markedly below the 
ductile brittle transition temperature because of major 
decrease in ductility. 

The effects of increasing temperature on erosion 
have been investigated in more details during the last few 
years . Depending on the material being eroded, the tempe- 
rature range, and the environment, different effects due 
to temperature have been reported by Young and Ruff [40]. 

The erosion rate has been found to increase with increase in 
temperature [41*42] . Rirtnie et al. as given in Ref „ [25] ' have _ 
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reported increased erosion for aluminium and stainless steel 
at higher temperatures. This effect, however, had a compli- 
cated dependence on the impingement angle. Decrease in 
erosion rate with increasing temperature has been reported 
for certain alloys also [44]. 

levy and Chik [45] have experimentally shown that 
erosion rate depends on the hardness of the impacting parti- 
cles. They have reported an increase in erosion rate upto 
a Viewer's hardness of about 700 kgf/mm 2 and thereafter it 
remains more or less constant, for both normal as well as 
o bl ique impact s . 

1.3.2 Brittle Erosion 

Because of lack of ductility, brittle materials are 
highly susceptible to micro-fracture or localised surface 
cracking during solid particle impact. Under continual 
multiple impacts, microfractures formed in brittle solids 
result in appreciable material loss by erosive wear. The 
formation of micro-cracks in brittle solids under such 
situations is a complex process and depends on the physical 
properties of the target and the projectile materials and 
the projectile velocity. Attempts to describe the erosive 
wear behaviour of brittle materials have tended to be 
semi-empirical. Several attempts have been made to use 



quasi-static indentation results for evaluating the extent 
of damage produced dxiring dynamic particle impact [46,48, 
50,56,63]. This approach has been adopted by Sheldon and 
Finnie [46] to analyse the erosive wear phenomenon of 
brittle materials. They assumed that the interactions 
between the impacting particle and the target are perfectly 
elastic a.nd Hertzian contact stresses are produced in the 
target during the course of indentation under the impact. 
These stresses result in the formation of Hertzian cracks 
which initiate from pre-existing flaws that lie just out- 
side the area of contact between the particle and the 
target surface. The erosion loss or material removal 
takes place because of the interaction of these cracks. 

The load at which the crack propogation occurs is related 
to the distribution of surface flaws through Weibull 
statistics. Using momentum balance, the maximum depth 
of indentation z^ was obtained. The final cracked space 
is delineated by radius of curvature r f which is ass um ed 
to be equal to that of the effective local radius of 
curvature of the indenting ball producing cracking when 
the ball penetrates to z . Through these the maximum 
fracture radius a f (the fictitious contact radius of a 
ball of radius r f ) can be obtained using the approximate 
geometrical relation 
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The volume removed is taken approximately equal to 
the volume of the spherical cap defined hy r (or a^) and 
z . Thus, the final expression of eroded volume V in terms 
of impacting particle radius r^, the impact velocity u 
and ¥eihull constants m and o Q is expressed as 


V 


c 4 



n. 


u. 


pxm 


(1.14) 


where the exponent n-^ for spherical particle is given hy 


= 3(m - 0.67)/(m-2) , 

and for angular particles is given by 

n 1 = 3.6 (m - 0.67)/(n-2) . 

The exponent n^ for both spherical as well as angular 
particles is given by 


n^ = 2*4(m - 0.67)/(m-2) . 


For particles much stiffer than the target, the constant 
c^ is given by 


■R. 


0 . 8 ( m+1 ) 1.2(m-0.67) 


'4 


jsl 


-3 


l/(m-2) 


(1.15) 
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■where B^ is the modulus of elasticity of the target material 
and Pp is the density of the particle material. Theoretically 
predicted exponents n^ and n 9 have been found to he in satis- 
factory agreement with the experimentally obtained values 
for several brittle materials such as glass, MgO, graphite 
etc. [46]. Sheldon- [47] compared the experimental and 
theoretical values of c^ and again found reasonable agreement 
between the theory and experiment. The agreement, however, 
was not as good as that for the exponents n^ and np. 

Inherently, this elastic-brittle treatment leads to 
the conclusion that erosive wear of brittle material depends 
on the fracture toughness [4-8] and on the size and the 
distribution of pre-existing flaws [46,49]. Investigations 
on the behaviour of brittle solids when quasi-st at ically 
loaded with small enough sphex'ical indent ers have revealed 
that the plastic flow does occur in brittle solids [48-53]* 

Oh et al. [49] have shown that transition from fully elastic 
crack formation to elastic-plastic indentation occurs, even 
in very brittle materials such as glass, when the diameter 
of the indenting particle decreases below a critical value 

d . The morphology of fracture formed in brittle materials 

IP ^ 

during impact can be divided into two classes, depending on 
whether the impacting particle is ’blunt’ or ’sharp' [48,54]. 
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Blunt particles, typified -by spheres, result in the 
formation of cotie shaped cracks (Hertzian cracks) -which 
initiate from pre-existing surface flaws that lie in the 
tensile stress field outside the contact circle between 
the particle and the target surface [54,55]* The entire 
nrocess of Hertzian crack formation is an elastic one that 
can be described by linear fracture mechanics. Bor quasi- 
static impart, lawn and Wilshaw [48] and Frank and lawn [56] 
have developed the requisite fracture mechanics equations 
for crack growth, langiton and lawn !"57] ha ve shown that 

s' 

the size of the flaw that initiates Hertzian cranks plays 
a secondary role in the crack growth. 

Because Hertzian cracks initiate from pre-existing 
surface flaws, a minimum threshold stress (correspondingly 
a minimum impact velocity) is required for its formation 
[58,59]. The threshold stress foir crack formation is 
governed by the diameter of the impacting particle and the 
critical stress intensity factor of the target material. 
Vie&erhorn and lawn [59] have obtained a theoretical expre- 
ssion showing the offoct of particle velocity, particle 
material density and particle size on the strength degra- 
dation of brittle material such as glass when impacted by 
gteol and tungsten carbide spheres. The critical load for 
crack growth is determined from the size of the critical 
flaw present on the target surface. 
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Knight et al. [60] have studied the crack systems 
produced due to impact of steel spheres (0.8 to 1.0 mm 
diameter) and tungsten carbide balls (~C.4 mm diameter) 
at different velocities on pyrex and soda-lime glass. 

They found that the crack systems produced dtie to impact 
of particles are similar to that obtained during static 
indentation. They further observed that a large number 
of fine splinter cracks are initiated beneath the impact 
site and propogate into the surface perpendicular to the 
maximum slip stress trajectories. 

The sharp particle indentation (with cone, pyramid 
etc. ) is characterised by a plastic flow at the contact 
zone which results in the formation of two types of crack 
systems ~ radial and lateral [50,54, 61-68]. The first 
type is a radial crack system oriented primarily perpendi- 
cular to the target surface, while the second type is a 
lateral crack system oriented parallel to the target 
surface [62]. The system of radial cracks is produced 
during the loading part of the impact cycle, whereas the 
system of lateral cracks is produced as the particle leaves 
the target surface [50,62]. Using e3.astic -plastic analysis, 
lawn and Evans [69] have obtained a theoretical estimate of 
the force for radial crack formation, further, they have 
shown that the critical force as well as the critical crack 
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size, for producing radial crack, depend on the hardness 
and the critical stress intensity factor of the target 
material. Tho threshold velocity, to produce the critical 
force which initiates the crack formation during sharp parti- 
cle impact [66], is considerably lower than that for the 
olunt particle impact. 

Swans et al. [63] have assumed that a sharp sphe- 
rical particle penetrates into the target surface without 
distortion and the dynamic pressure is setup when the 
particle first hits the surface. This gives rise to a 
contact stress field. Using one-dimensional impact analogy, 
they obtained the expression for the size of the radial 
crack formed during impact in terms of velocity, density, 
size of the impacting particle and the critical stress 
intensity' 1 factor of the target material, from energy 
balance considerations Wiederhorn and lawn [66] have 
developed a relation between the contact force, hardness 
of the target material and the maximum depth of penetration 
during impact, further, they have also expressed the size 
of the radial crack formed during impact as a function of 
the target hardness, the critical stress intensity factor 
of the target material, the mass and velocity of the impact-, 
ing particle. Tho final expression for tho radial crack 
size obtained by Evans ot al. [63] is quite different from 
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that obtained by Wiederhorn and lawn [66]. However* both 
theories suggest that crack formation depends on critical 
stress intensity factor of the target material. The theory 
of Evans et al. [63] suggests a stronger dependence on the 
velocity of impacting particle. 

The plastic zone beneath the contact point as well 

as the lateral crack depth during sharp particle indent a-' 

tion are proportional to the impression diameter [50]. 

Assuming that the depth at which lateral cracks are formed 

is proportional to the plastic zone depth* Evans et al. [63] 

have obtained the relationship between the depth of lateral 

crack z n , particle size d , -particle velocity u . * particle 
1 p pun * 

density and target hardness EL as 


( Z l)2 

v d ' 

P 


u 


-pirn 


r ~t 


( 1 . 16 ) 


The radial cracks are primarily responsible for 
strength degradation, whereas the lateral cracks are 
responsible for erosive wear of brittle materials [50,62], 
The erosion rate of brittle materials in elastic -plastic 
damage regime is sensitive to both the size distribution 
as well as the velocity of the impacting particles [70]. 
Experimental studies of Scattergood and Routbort [71] have 
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revealed that the velocity exponent in the erosion rate 
equation depends on the size of the impacting particles. 

The velocity exponent decreases with increase in the size 
of the particle. They found the velocity exponent to 
decrease from 3.4 to 2.6 when the particle size was incre- 
ased from 23 to 270 pm. 

Tabakoff et al. [72] have revealed that the angle of 
impingement corresponding to maximum erosion depends not 
only on the target material but also on the impacting 
particle material. 'The erosion damage is more dependent 
on relative target - projectile properties than on the size 
and imping-ement velocity of the particle for a given target 
projectile particle combination [73]* 

Mehrotra et al. [74] have considered three mechanisms 
of material removal under impacting particle : material 
chip-off when cone cracks form on the surface; material 
removal due to interaction of cone cracks on the flat 
surface and material removal duo to chipping of exposed 
cone frustums as well as the underlying rough surface. 

The proposed erosion model was verified at four different 
particle velocities and good agreement between experimental 
and theoretical prediction was obtained. 



Wiederhom and Hockey [75] have applied a dimen- 
sional. analysis of the variables that have significant 
effect on erosion. A multivariable linear regression 
analysis was used to fit the data to the dimensional 
anal^/sis. They have shown that the dependence of erosion 
rate on hardness and critical stress intensity factor is 
greater than predicted by erosion theories [63,66]. 

1.4 Present Work 

The literature survey reveals that very little is 
known about the abrasive jet machining process. Most of 
the relevant papers seem to deal with the problem of 
impact erosion. In the present work an attempt has been 
made to carryout theoretical as well as experimental study 
of the abrasive jet machining process. 

Material removal in AJM occurs due to erosion caused 
by the Impact of high velocity abrasive particles carried 
in a fluid stream onto the workpiece surface. The erosion 
studs?" in such situations involves the determination of 
number, direction and velocity of abrasive particles strik- 
ing the work surface from fluid flow conditions and then 
the estimation of material removed from the workpiece. In 
the first phase of the study the mechanics of air-abrasive 
mixture flow through convergent -straight nozzle as well as 
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jet flow have been developed. From two-^phase flow analysis 
a theoretical, model has been proposed to estimate the velo- 
city of abrasive particle in the air-abrasive jet. 

During solid particle impacts lateral cracks are 
produced in brittle materials below the impact site as a 
result of residual stress field. Erosion loss in brittle 
materials occurs due to intersection of lateral cracks from 
adjacent impacts. On the basis of quasi -static indentation 
approach and using semi-empirical fracture mechanics relation- 
ships, the length of lateral crack and the depth at which 
lateral, cracks are formed have been expressed in terms of 
impact parameters. Subsequently, in the second phase of 
the study, a theoretical model for material removal rate 
has been developed with the informations of lateral crack 
sizes . 

In polishing, deburring and other finishing opera- 
tions the weight /volume of the eroded material is of 
importance, while In micro-drilling and cutting operations 
it is the erosion depth (depth of penetration) which is 
more relevant. The third phase of study involves the experi- 
mental investigation of different parameters affecting the 
erosion volume, erosion depth and eroded cavity diameter etc. 
Eor this an experimental Abrasive J et Machine has been 
designed and fabricated. The setup provides precise 
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control of the input parameters such as air and abrasive 
flow rates, nozzle pressure, mixture ratio and stand-off 
distance. 

In AJM the nozzle, through which the air-abrasive 
mixture is expanding, is one of the very important compo- 
nents. The effectiveness and performance of the Abrasive 
Jet Machine depends on the geometry of the nozzle. The 
high velocity air--abrasive mixture flow through nozzle 
causes wear or erosion of the inner surface of the nozzle 
which results in changes in the geometry of the nozzle. 
Eozzle wear is therefore of interest not only from the 
point of view of cost but also from the point of view or 
characteristic performance of the unit. In the fourth 
phase of the study an attempt has been made to carryout 
experimental study of nozzle wear. The parameoeiS that 
have been considered are nozzle length, nozzle diameter, 
entrance angle, mixture ratio and particle size. 

In the last phase of the work, the theoretical as 
well as experimental results are discussed, in detail in 
order to understand the mechanics involved in the process. 
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Table 1.1 : Typical Applications of Some Commercial 
Abrasives [After Bhattacharyya( 4) ] 


Abras ives 

Grain Sizes 

Machining Operation 

Aluminium Oxid e 

( Al^O ) 

2 3 

12,20,50 pm 

Cutting and grooving 

Silicon Carbide 
(SiC) 

25,40 jim 

Cutting and grooving 

Sodium Bi- 
carbonate 

27 pm 

light finishing opera- 
tion below 50°C 

Dolomite 

76 pm 

Etching and polishing 

Glass Beads 

0.655 to 

1 . 27 mm 

light polishing and 
fine deburring 


Table 1.2 ; 

Dimensions and 

of Nozzles and 

charyya(4) ] . 

Nozzle life for Various Types 

Nozzle Materials [After Bhatta- 

Nozzle 

Round Shape 

Rectangular Shape life of 

Mat erial 

Nozzle Dia- 

Slot Dimension, mm Nozz3_e, 


meter, mm 

Hours 


Tungsten 

Carbide 


(WC ) 


0.2 to 0.1 0.075 x 0.5 to 12 to 30 

0.15 x 2.5 


Sapphire 


0.2 to 0.8 


300 
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rate for different work materials 
io) cemented carbide end ib) glass. 
Abrasive : aluminium oxide jparticle size: 
40 pm f After Sorkar and Pandey(5)J . 
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F in 1.9 Volume remove! cs o function of the hardness 
y f nr rn eta Is eroded at 9=20 and velocities of 

76 m /sec and i37 m/sec. d p : 250 pm Si C 
[After Finns© ef ol. (10)] 



CHAPTER 2 


THEORETICAL ANALYSIS OF ABRASIVE PARTICLE VELOCITY 
2 « 1 Int roc 1 net i on 

Eor many years engineers and scientists have been 
interested in gas -solid particle flows which arise in many 
industrial applications . To mention a few gas-solid 
particle flow phenomena are important in sedimentation, 
pipe flows [76-78], fluidized beds, transport processes 
r 79,80], propulsion and combustion [81-83], abrasive clean- 
ing devices (sand blasting units), abrasive jet machining 
etc. In recent years, gas-solid particle mixture flow 
through rocket nozzi.es has received much attention [84-89]. 

Gas -solid particle mixture flow is a particular 
class of two-phase flow in which small solid particles are 
suspended in a gas which acts as a carrier phase. The 
continuum mechanics approach similar to that for single-phas 
flow has also been found most powerful and suitable method 
to deal with two-phase flow system [90]* The basic diffe- 
rence between single-phase flow and particulate-phase (gas- 
solid particle) flow lies in the mechanism of information 
transfer between the elements of the particulate-phase. 

There is no analog for pressure in a particle cloud and no 
information transfer exists between particles except through 
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the conveying phase. Moment um and energy interchange occur: 
between two phases mainly by viscous interaction and heat 
transfer, other sources of interchange, like mass transfer 
by evaporation and condensation or chemical reactions have 
not been considered. 

In this chapter an attempt has been made to bringout 
the basic characteristic features of gas-solid particle flow 
using certain simplifying assumptions. Various concepts for 
understanding gas-solid particle flow through a convergent- 
straight nozzle and estimation of particle velocity inside 
the nozzle as well as in the jet region (i.e. after exit 
from the nozzle) have been presented. 

2.2 Basic Assumptions 

for the treatment of particulate flow using fluid 
dynamical equations, the particle phase is assumed ao 
continuum. All particles are assumed to be uniform sphe- 
rical particles of identical physical properties. The 
particles are assumed to be sufficiently large so that 
Brownian motion of particle can be neglected. It will be 
further assumed that the particles are sparsely distribyued 
so that flow fields around the particles are not affected 
by other particles^ Thus, the particle cloud may be assumed 
to be pseudofluid and fundamental equations of motion can 
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also be applied for the dense phase. Also, the solid 
particles are assumed to be uniformly distributed at each 
cross-section although it is understood that they are in 
suspension due to turbulence and interaction exists between 
components of the mixture. 'The motion is assumed to be 
one-dimensional and steady and the effect of turbulence has 
been considered only through characteristic parameters. 

Mass or energy loss from the system is neglected, i.e., 
there is no mass exchange between the phases. The fluid 
phase is assumed to behave as a perfect gas and that it is 
inviscid except for its interaction with the particles. 

The body forces, viz. gravitational, electromagnetic, 
buoyancy etc. are neglected. 

2.3 Aerodynamic Response Time 

The density of the solid particle material is much 
greater than the gas density. Thus, solid particles can- 
not accommodate the rapid change in velocity (acceleration 
or deceleration) and tend to slip through the gas as they 
accelerate or decelerate. This particle slip velocity gives 
rise to a drag force acting on the carrier fluid and an 
equal and opposite drag force acting on the particle phase. 
Following the particle path and neglecting Basset force, 
the equation of motion for spherical solid particle of dia- 
meter d and material density p in a uniform stream of gas 

r P 

Of density P and velocity u is given by [90-92] 

§ & 
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f a ! ? p av = c d [ * ? s ( vV’ V + ? a p 5 s 3t^ 


+ * tf d p p ”g af~ <VV 3 ’ C2 - 1) 

P 

where A is the projected area of the solid particle 
P 

perpendicular to the flow, is the corrected drag 
coefficient, u^ is the velocity of the solid particle 
and tp is the time for the motion following particle path. 

The first term on the right hand side of equation 
(2.1) is the viscous drag terra. The second term is due to 
pressure gradient in the fluid surrounding the particle. 

The third term is the force to accelerate the apparent mass 
of the particle relative to the fluid. 

Effects of pressure gradient, apparent mass, and 
Basset force become important only when the density of the 
fluid is comparable or higher than that of the solid parti- 
cles [90]. In the present case, the ratio of densities of 
the solid particle to that of the fluid is of the order of 
10 3 . Thus, the terms due to pressure gradient and apparent 
mass are also insignificant and equation (2.1) reduces to 



( 2 . 2 ) 
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The projected area Ap for spherical particle 
■perpendicular to the direction of flow is 




T d 2 
4 p 


(2.3) 


Stoles [93] has given the drag coefficient C^ g 
for spherical particle moving in fluid as 


C 


24 


ds 


Iff. 


Re 


(2.4) 


Here the Reynold's number, is given by 


E 


Re 




(2.5) 


where "u is 

8 


the viscosity of the fluid 


phase material. 


Stokes drag law holds good for < 1. For flows 
at higher Reynold's number, Gilbert et al. [94] has intro- 
duced a correction factor f in the drag equation given by 
Stokes [93] and expressed the corrected drag coefficient 


as 


J a 


= f C, 
c ds 


where the correction factor f was given as 

o 


( 2 . 6 ) 
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f 

c 


, -0.85 

(28 E- Re + 0.48) 



(2.7) 


Using equations (2.4), (2.5) and (2.6), the corrected 
drag coefficient 0^ can he expressed as 


C 


d 


24 h f c 

s e 4 P <vV 


( 2 . 8 ) 


Equation (2.2), after substitution of equation (2.3) 
and (2.8), becomes 
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( 2 . 10 ) 


which is known as the 'time constant' for momentum transfer 
due to drag force, equation (2.9) becomes 
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equation (2.1l) can be written as 


^• u p ^p 

% + i 


u g 


T 


(2.15) 


p 


Por t — o and u = o, solution of equation (2.13) yields 
P ' P 


u s = U g (1 


tp 

e P ) . 


(2.14) 


For t 


P 
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u 


= u_ 


Lg-iJ. 

e 


(2.15) 


This time 1 is defined as the particle relaxation 
1? 

"h j_t n 0 or 3 , 6 Tocl jrD-stniic j^osponsB xiuiib# X"t is iliG iiiiiB 3?6(pu.ii*Bcl 
for a particle released from rest in flowing stream to 
achieve 63 percent of the mean free stream velocity T 90 9 95 ] - 
t:_P is the characteristic time of the flow system and can 
be expressed as 


where 1 is the characteristic length of the flow system 
and U is the characteristic velocity. Stokes number 

S_ r is defined as the ratio of aerodynamic response time 

J i 
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to the characteristic time of the flow system. Thus, 

. (2.17) 

xv T ^ 

When Stokes numher is very small, the particles 

have ample time to respond to property changes in the 

conveying phase, i.e. the particles adjust to the system 

very quickly and their characteristics are governed by 

local conditions. This state of flow is described as 

'equilibrium flow'. On the other hand, when the Stokes 

number is very large, the particles are not able to respond 

to the changes in gas flow. Such a flow is generally called 

‘frozen flow' and the particle 'Characteristics are governed 
by their initial history. 

2.4 Behaviour of Particulate Plow 

The relative volume and mass of a component in a 

gas-solid particle mixture are quantified by the volume 

fraction and bulk density. The volume fraction is the 

volume of a phase per unit volume of mixture while, the 

bulk density is the mass of a phase in Linit volume of mixture. 

Por gas-solid particle mixture, if Q denotes the volume 

fraction of solid particles then the bulk densities of 

particle phase p and fluid phase can be expressed as 

P s 



p 

p 


(2.18) 
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and 

D. = (1 - Q } P, . (2-19) 

& & 


where P 

P 

and fluid 


and P „ 
g 

phases 


are the material densities of particle 

respectively. When Q is small, p ~ P . 

§ & 


'The mixture rat io a as indicated earlier is defined 

as the ratio of the mass flow rate of the solid particles 

« 

Np to the mass flow rate of gas and solid particles mixture 

» 

M. Thus, mixture ratio can he expressed mathematically as 


a 



Alt ernat ively , 


and 


Tip = a M , 

M = (l - a) M . 

g 


( 2 . 20 ) 


( 2 . 21 ) 


( 2 . 22 ) 


The coefficient of viscosity for the mixture of gas 
and solid particles is not the same as that of the gas phase 
alone. It is also not a constant because the density of the 
particle phase changes TPO]. Sproull [96] has measured the 
viscosity of dusty gases (air and dust particles of lime- 
stone and talc) using rotating viscometer and found the 
viscosity of mixture to he less than that of clean air. 
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There appears to be some controversy about the viscosity 
of gaS'-solid particle suspension [90, 96-99]- Although 
there is disagreement as to whether the viscosity increases 
or decreases when solid particles are added to the gas, the 
viscosity of the mixture jj. » according to Einstein [See 
Reference (90)], can be expressed as 

= \i g (1 + 2.5 Q ) • (2.23) 
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the viscosity of the gas phase alone. 
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(2.25) 


where u is the viscosity of the particle phase and u is 

1 p " - * 1 g 

the viscosity of the gas phase in the mixture. 


Eor very dilute suspension, 


h « p g ’ (2 - 26) 

and 

/v o , (2.27) 
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i.e., viscosit y of the gas phase is approximately equal 
to the viscosity of gas phase material and the particle 
phase viscosity may he neglected. 

2.4.1 Squat i on of St at e 

fhe equation of state for the gas phase is 

p g = R e p e f e " V^’sV (2 ' 28) 

where P , is the partial pressure due to gas phase, R is 

g S 

the gas constant and T is the temperature of gas phase. 

§ 

The total pressure P from the perfect gas law is 

P = B. P' T , (2.29) 

g g g 

and 

P _ = (1 - Q) P , (2.30) 

O 

where P = (P + P ) . 

e> y 

H enc e , 

P r , = Q P = • (2.31) 

P S S e> 

This indicates that when Q is small the partial pressure 

P_ of the particle phase is negligible. This fact is also 
P 

supported by saying that the particle phase has no analog 
of pressure. The equation of state for the mixture can, 
therefore, be expressed as 
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g g 


( 2 . 32 ) 


2.4.2 C ont limit y Bquat ion 

The continuity equations for the case of two-phase 
(gas-solid particle) flow are similar to those for a single- 
phase flow with the exception of source terms. These equa- 
tions for uniform particle sizes are [90] 


3 p p 3 Pp u p 

at" + 'ax?” 



1=1, 2, 3 


(2.33) 


for the particle phase, where t is the time, x^ is the 
ith component of the space coordinate, u^ ± is the ith 
component of the velocity of the particle phase and 


+ ^gi 
at a x i 
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p 


(2.34) 


for the fluid phase. Here u . is the ith component of the 
velocity of fluid phase. The term Y normally called as 
source function of particles, represents the mass addition 
to the fluid phase per unit volume of mixture caused by 
phenomenon such as evaporation, condensation etc. If there 
is no change of phase, as in the present case, the source 
funct ion ¥ = 0 . 
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2.4.3 Limit of 0 ont inuum Appr oximat ion 


When solid particles are introduced into a fluid 
stream there is a velocity distribution around the particle 
due to the relative motion. When the relative speed is 
high enough, the existence of wake is expected. Therefore, 
the concept of continuum mechanics when applied to the 
case of particulate suspension requires examination of the 
nature of fluid flow around the particles. 

In the absence of shear motion, body force and chemical 
reaction, the equation for conservation of momentum of 
suspension of identical solid particles can be written as 


du 


J g at 


jB3 


du. 


’p dt 


lEi 


D 


8 L 

3 x. 


(2.-35) 


where t and t are the time for the motion following fluid 
o P 

path and solid particle path respectively. In this equa- 
tion the particle phase is treated as a true continuum, 
that is, not only the gas imparts its momentum to the 
particles but the cloud of solid particles can also transfer 
its momentum to the gas. This is obviously not always the 
case. Even in laminar flow range of relative motion, Fross- 
ling [100] have shown that the boundary layer thickness 6 
for flow over a sphere can be evaluated from 
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-S 5 . (2.36) 

P P g 

Therefore, if cloud of particles is to exert sufficient 
drag on the gas, the interparticle spacing must "be smaller 
than (26 + d^) . 

The number density of solid particle n is 

Jr 

defined as number of solid particles per unit volume of 
mixture. A cube of sides (d +26) will contain eight 
particles one at each corners with one eighth contribution 
to the considered cubical volume of mixture. Tims, the 
effective number of particles per unit volume of mixture 
is one. Hence, the number density of solid particle is 




26 + d n< 26 

P ' 


(2.37) 


This relation loads to 



( 2 . 38 ) 


for successful transfer of momentum, of mass motion between 
tbo two phases. Here the solid particles may also accelerate 


the fluid phase in a similar manner as gas accelerates solid 
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particles. ¥hen equation (2.38) is not satisfied, the 
fluid phase (gas) simply acts as a free stream dissipating 
the kinetic energy of the particles in their wakes. This 
physical state leads to the conclusion that outside the 
Stokes 1 law range, drag coefficient of solid particles 
depends on the concentration of clouds "besides the 
Reynold's number of relative motion. 

To account for all situations, Soo [90] has shown 
that equation (2.35) should he modified as 


du. 


du 


v 


dt 


Jai 


+ 


K m 


p at 



(2.39) 


with the 'effectiveness' K = 1 for the case when solid 

Li 

particles are accelerated "by the gas, and E ffl < 1 for the 
case when solid particles are decelerated by the gas. 

For a dense cloud of solid particles, when negligible 

dissipation in wake is encountered, E = 1. For a lean 

suspension, dissipation in wakes without increase in 

kinetic energy of mass motion of the gas gives E ^ o. 

Effectiveness F also denends on the concentration of the 

m 


stispension. 
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The exact reoresentat ion of K was idealized using 
- m 

a discontinuous representation as shown in Fig. 2.1, 
although a smooth transition (dotted line) is expected [90], 
It is further noted that the momentum transfer in particu- 
late phase depends not only on the drag coefficient but 
also on the effectiveness I. which in fact depends on the 
concentration of the suspension. Therefore, the time 
constant of momentum transfer B depends on the concentra- 
tion. of the suspension in addition to the dependence of 
drag coefficient. With this understanding of the continuum 
approximation, the basic equation of motion for two--phase 
flow system can be modified using E for the particulate 
phase. This modification accounts for the fact that slowing 
down of solid partic3.es does not necessarily give rise to 
increase in static pressure of the gas. 

2.4.4 Squat ion of Hot ion 

Fai [lOll has shown that under continuum conditions, 
particles in pseudofluid assume stress-strain relationship 
similar to that of a gas. The body forces in gas-solid 
particle flow may arise because of gravitational forces, 
electromagnetic forces etc. The overall momentum equation 
and the particle phase momentum equation as expressed by 
Soo [90] are 
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dUp-i du pi 

° + K p. - 


'g dt T m dt 


P 


9P , ^ 

5x ± p g Dgi 


+ [ n ( a ) . . + \i D 6..] 

3x. LK g g 01 Tg ? g oi- 
3 <L 


(2.40) 


and 


dUpi 


p dt 


P 


•° B(u . -u . ) + P f, . 
p gi pi' p Tapi 


+ 53T Wji + M S 2 B p 6 u ] - (2 ' 41) 


■where is the ith component of the "body force per unit 

mass of the fluid material, ^ is the ith component of 

the external force acting on spherical particle per unit 

mass of particle material in general fluid field besides 

the drag force and P is the total static pressure in the 

fluid. The terms p, , and \x represent the -viscosity of 

S>2 

the fluid and particle phases respectively including the 
bulk viscosity of the respective phases i.e., 

% S g " 3 ' 


M- 


P 


2 


5 


y ' 2 

h " 3 ^p 
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where ^ 
particle 


and £ are the hulk viscositjr of fluid and 
g P 

phases respectively. In equation (2.40), 


a 


dt 


3 o 

Ft + u gj 3 xT ? 


d_ 

dt 


P 


= -,rr + U . *r , 

dt po dx. ’ 


and the deformation tensor (a and dilation D fo3 

S a 

the fluid phase are 


^ A g^di 


5u 


gi 


dx. 

J 


5u ii 

az ± 


D 


1 (a ). . 

2 v gkk 


® u gk 

3x. 


In equation (2.41), the deformation tensor (A).. 

x Da 

and dilation for the particle phase are 


(a ) = + iii 

v p j i ox. ox. ’ 

x- <j 1 21 


3u- 


D 


(a ) , i 

2 p kk 


pk 


v k 


Here i, j, k = 1, 2, 3 ; Krone eke r delta 6.. = o for 

J 

i ^ 3 and 6.. = 1 for i = j. 

ci 
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2.4.5 Energy Equat i on 

The energy equation of a two-phase (gas-solid 
particle) flow system includes the relationship of 
exchange of energy between the fluid phase and the particle 
phase. Eor particulate flow, that is for both the fluid 
and the particle phases, the overall energy equation for 
mixture can be written as [90] 


_c L 
dt 


u 


gi 

2 


u. 


c T ) 
Pg S 


+ P 


P dt. 


Pi 

2 


c T ) 
PP P 


3P 

^t 


5 3 9 ^p 3 

+ HTj" k g sU + 3x7 k p 3x~ + 3x7 


+ ^ Vu ] + W Vji + h 2 Vji ] 


+ r g + r p . 


(2.42) 


where c , T and k are the specific heat, temperature and 
P 

thermal conductivity. Suffix g represents the fluid phase 
while p represents the particle phase. The term P. includes 
the heat source and the radiations per unit volume of the 
mixture and includes the heating of particles due to 
viscous dissipation. 
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2.5 Particle Velocity in the ITozzle 

In the present analysis, under the assumptions 
made in Sec. 2.2, one-dimensional steady state friction- 
less adiabatic flow through convergent -straight nozzle 
with dilute concentration of solid particles has been 
considered. With the introduction of mixture ratio a, 
the continuity relation (Sec. 2.4) gives the flow rate of 
gas phase as 


Pg u g = (1 - a) M = Constant , (2.43) 

and the flow rate of particle phase as 

Pp Up = a M = Constant , (2.44) 


where is the cross-sectional area of the flow passage 
of the nozzle. 


For steady state frictionless, adiabatic flow the 


terms 


0 


0T 


g 


01 , 


fx7 0zT 5111(1 It of equation (2.42) 
3 D 


ox^ g 0X.. » 


become zero. For inviscid fluid phase, neglecting heat 
losses (heat transfer due to convection, radiation through 
wall, heat transfer between particle and wall due to impact 
and heating of particles diie to viscous dissipation) the 
overall energy equation of the mixture for the gas-solid 
particle flow through the nozzle reduces to 
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u. 


p g u g dx ^ c pg T g + 2 + p p U p dx ^ C PP T P + 2 ^ ° 


(2*45) 

Differentiating -with, respect to x, and using equations 
(2.43) and (2.44)» equation (2.45) becomes 


du r 
'g dx 


dT_ _ dtu 

S , { K \ r„ P 


dl. 


U V"” + C 5 * + -) Tu 5 — “ + C 

- pg dx v l~cr L ^> dx — ' " J 


P- 

PP dX J 


= 0 . 


( 2 . 46 ) 


For steady state, invis c id , non- reactive and un- 
charged gas-solid particle flow equations (2.40) and (2.41) 
become 


du 


p u , 
g g dx 


g 


du, 

K m |0 p u p dx 


dP 

dx 


(2.47) 


and 


u 


P 


dUp 

dx~ 



(2.48) 


Equations (2.43) and (2.44) are the continuity 
equations, equation (2.46) represents the energy equation, 
equations (2.47) and (2.48) represent the momentum equa- 
tion, for the gas-solid particle mixture flow through the 
nozzle . 
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Rearranging, equation (2.46) yields 


dT^ 

dx" 


, u g > ^ 
( “ } dx 
Pg 


a 

l~a 


u d 

r(-E-) ~ 

u c ' dx 
pg 


^ r5> + ( 2ev, Ez, 

+ { c dx J 
Pg 


(2.49) 


Substitution for p from equation (2.43) into 

o 

equation (2.32) yields 


P . M iVi 


(2.50) 


Differentiation of equation (2.50) with respect to 
gives 


dJP 

dx 


( 1-a ) M R [(*--—-) n 
g L u A dx 
g n 


dT, 


T du T 

( — £_) —S _ 

v , 2 ' dx 1 . 2 


Vg 


A u 
n g 


) — ]• 
‘ dx J 


(2.51) 

Substitution of equation (2.43) into equation (2.51) and 
after simplification yields 


dp r dT s d ^ e ^ dA 

ta = p g R g [ if - ^ dir) ■ ( 2 - 5a > 


Using equations (2.43) and (2.44), equation (2.47) can 
be written as 
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dP 

dx 


Solution 


cl Up, „ dUv, 

= - p u + Z (-/-) v-£] . 

g g L d.x m 1-a dx J 


of equations (2.52) and (2.53) gives 


(2.53) 



du r 

' n ' 


u. 


dx 


I, dA. 
+ (-£) — -Si. 
dx 


u. 


(- 6 ) 


du. 


Vi 


a 


) ^ 

<x R. 


du. 


■D 


(2.54) 

dT g 

Elimination of from equations (2.49) and (2.54) and 

P 

substituting for R = pr-m gives 

8 P g g 


du 


g 




Vs* A u 


_ JL_) 

c ' 


duj 

dx 


(-- 

u 


u. 


P c 


yv_ 

'm 

T 


dT,, Tg dA 
; dx * dx J 


g 


PS 


(2.55) 


Rearranging , equat ion (2.48) gives 


clx 



1) 


(2.56) 


2.5.1 Equation of Heat Transfer for Flow Over 
Spherical Particle 


G-ehhart [102] has given an empirical equation for the 
Husselt number correlating data for Reynold’s number 
upto 70,000 and Prandtl number Fp r upto 400 for the case of 
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heat convection from the surface of a spherical particle 
immersed, in a fluid stream* This equation is 



_ 0.5 1/3 

2 + 0.6L IL 
Re Pr 


(2.57) 


From the definition of the Nusselt number [102] equation 
(2.57) can he rearranged as 


h 


g Vu 


( 2 + 0.6 1 

d' 


0.5 

Re 


1/3 
H lr ) 


(2.58) 


where h is the coefficient of heat transfer between the 
particles and the gas and k is the thermal conductivity 
of the gas phase. 

The heat exchange between particle and gas during 
time dt is given by 


h A sp w at 


“p °pp (_ d V - <2 - 59) 


where dt is the time taken by the particle to travel a 

distance ds and A is the surface area of the particle. 

sp 


Thus, particle velocity u^ 


is 


= u 


ds 

dt 


( 2 . 60 ) 



f 


Mass m and surface area A of the particle 


•o 


sp 


are 


m 


P 



( 2 . 6 : 


and 



( 2 , 62 ; 


Dividing equation (2.59) by ds and substituting 
for h, 4r , m and A„^ from equations (2.58), (2.60), 

CL 0 p Sp 

(2.6l) and (2.62) gives 


T -T 

dT„ 6 k s g p n - - , 

™JL — fp 1 a 6 t\t® * 5 i\rl/3 \ 

dx - p a 1 ®e h?r > ' 

P P PP 1 


(2.65) 


2.5.2 Dimensionless Variables 


It is more convenient to use the gas -solid particle 
flow governing equations in dimensionless form. The common 

method of non-dimensionalizing the nozzle flow governing equa/ti< 
in terms of Mach number is not suitable for, gas-solid flows-. 

For a two-phase (gas-solid particle) system this technique 
may give misleading results because the simple relationship 
for the velocity of sound in clean gas is no longer valid 
for a gas-solid suspension. Although, the basic equations 
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may still Toe non-dimensionalized in terms of the gas phase 

Mach number, care must he exercised in interpreting the 

results. It may he mentioned that the Mach number will he 

incoirect except for very low mixture ratio. To account 

for general situations of dense as well as dilute suspension 

the different variables of the flow governing equations are 

non-dimensionalized below with respect to the inlet flow 

parameters of the gas phase and the outlet diameter d of 

no 

the nozzle. The dimensionless variables are indicated with 
superscript * . 


* 

u 

S 

U g 

U J 

SI 

* 

“'81 

* 

P 

g 

0 

_ 

~ p * 

81 


rn* 
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T„ 
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- — , 

SI 

u* 

p 

_ 

T * 
gl 

•vf 

P 

P 

P P * 

" P ' » 
gl 

> (2.64) 

* 

P 

P 

= P x ' 

* 

" K'x 

* 

■v 

X 

~ d » 

no J 



where A „ is the area of nozzle cross-section at the inlet 
ri-L 

u P~_» are velocity, density, temperature 

g 1 bl b 1 ± 

and pressure respectively of the gas phase at nozzle inlet. 

# 

Substitution for M from equation (2.21) into equation 
(2.50) gives 
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.Reducing equation (2.65) to inlet condition gives 


u 


gl 


(1-cQ Hp 
K P 1 


( 2 . 66 ) 


Using equation (2.43) for inlet condition, and 
substituting for M fror. equation (2.21), gives 


P 


gl 


(1 - a) lip 


a u 

gl 



(2.67) 


2.5.3 Flow Governing Equations in Dimensionless Form 
for the Gas-Solid Particle Flow Through Nozzle 

Using equation (2.43) together with equation (2.64), 
the gas phase continuity equation in dimensionless form can 
be expressed as 

P* vl* 4 = 1 . ( 2 . 68 ) 

Through equations (2.44) and (2.64) the continuity equation 
for the particle-phase in dimensionless form can be written 
as 




a 

( l'-re) 


(2.69) 
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Similarly, equations (2.32), (2.49) and (2.52) in dimension- 
less form "become 


P* = P * I* 


(2.70) 


k c T , 
Pg gl 




du* du^ 

I* + u* —~1 


5 dx* 


e , * 
° dx 


C d'T 

/ c \ °PP ai P 

l-g c , * ’ 

pg dx 


(2.71) 


p; [ 

s dx* 


Tj du* 

(. J) 

u* dx* 


T* dA* 

(-£) -JS 1 
< ds * J 


(2.72) 


Substitution of P* from equation (2.68) into equation (2*72) 
gives 


/ A * 
£ n 


_ 

rn ** 

(if-) 

2i _ 

m * 

-j 

-dx* 

U 

g 

dx* 

< 

dx* 


Equation (2.55) rewritten in the dimensionless form and 
substitution of equation (2.70) yields 


rc*)[s(%a Ku *_ ( _i|i 


< du* c~ dl* T* dA* 

,( „) u * } — . k . — |]-(~§) — S] 

pggl P dx* pg dx* dx* 

-.2 2 ^ 


^ U 

[ 2jS + u * 

L + U g Vj 


“ « c t>g T «. ^ ax* yt'ax* 


r 2 p 

si £a • 


pg gi 


(2.74) 
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Dimensionless form of equation (2.48) becomes 


-O* ! JL 

■°1 1 u*’ 

p 


(2.75) 


Here B* is a dimensionless constant given as 


_ li_ d 
3. & no 

4 -,2 ~ 

d p u , 
P gl 


®Re < 28 ®Re 


0.85 


+ 0.48) . (2.76) 


Substitution of dimensionless variables from equa- 
tion (2.64) into equation (2.5) and using equation (2.68), 
the particle Reynold’s number can be expressed as 


( is Dai?. ) a . %) 


(2.77) 


The dimensionless forms of the cross-sectional area 

dA n 

At of the nozzle passage and its derivative — ~ for the 
n dx 

considered geometry of the nozzle may be expressed as 


2x *t an8 


1 \ 2 


) for z < f , 


* 4 o tanp 9 2x*tan(3 9 9 

< - [(1 - *=£ + 'atd -?r^ ] 


_ _ * * 

for 1-, < x < 1 9 , 
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(• ^) 2 
d nl 


for x* > i: 


(2.78) 


4tanj3- 


( 1 - " 


2x*tanP-, 


) for x* < 1- 


«“f 2 r „ d n2^ tan ^2 , & n2, 2x 

a* C( eC.7 tanp i + iL, " a*, ' 


for 1* < x* < 1* , 


0 for x* > 1* 


(2.79) 


where 


- < 2 ) 

2t an(3 


(2.80; 


d nl (ta J p 1 - tai^) d^ 

2 tan.p 1 tao-P^ tanP 2 ' 2 ' 2tan(3 2 * 

(2.81) 


(2.82) 


(2.83) 


as derived in Appendix I. Here d ^ is the inlet diameter, 
d n2 ds "^e diameter at the end of first convergent region, 
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d*-, is t’.,e dimensionless inlet diameter and d* is tiie 
nl n2 

dimensionless diameter at the end of first convergent 
region of the nozzle. The terms 1-^ and ±2 represent the 
dimensionless distances at the ends of first and second 
convergent regions measured from the nozzle inlet. The 
terms and ^ are semi-entrance angles of the first 

and second convergent regions of the nozzle. 

Equation (2.63) rewritten in dimensionless form 

"becomes 


dT 
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2.6 Particle Velocity in the Jet 


(2.84) 


The gas-solid particle mixture after expansion 
through the nozzle appears in the form of a jet at the 
nozzle exit. Because of the inertia effect, the solid 
particles do not attain the same velocity as that of the 
gas-phase at the nozzle exit. Thus, there exist a large 
slip velocity. It is possible to investigate the two-phase 
(gas-solid particle) jet flow in a manner analogous to the 
single-phase gaseous jet because of the large volume fraction 
of the gas~plia.se even when the initial concentration of the 
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solid particles is large. Por a mixture ratio a = 0.4, 
the volume ratio of air and solid particles is of the 
order of 10 3 « As the jet moves away from the nozzle, the 
volume ratio increases further due to the entrainment of 
surrounding air . 

Because of the friction developed on its periphery, 
the emerging jet carries with it some of the surrounding 
fluid which was originally at rest. This friction also 
causes the jet to spread outwards in the downstream 
direction and a jet Boundary occurs between the two streams 
which move at different speeds in the same general direc- 
tion. Such a surface of discontinuity in the velocity of 
flow is unstable and gives rise to a zone of turbulent 
mixing downstream of the point where the two streams first 
meet. The width of this mixing region increases in the 
downstream direction. This turbulence also causes mixing 
of the emerging jet with the surrounding fluid. Particles 
of the fluid from the surroundings are carried by the jet 
so that mass flow increases in the downstream direction. 
Concurrently the jet spreads out and its velocity decreases, 
but the total momentum remains constant. 

The theoretical as well as experimental investiga- 
tions for single phase jet [103,104] have revealed that the 
velocity profile at any section of the jet is of universal 
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nature. Because of the mixing of the surrounding fluid 
the jet velocity decreases. Thus, there are three distinct 
regions in a single-phase free jet; (i) potential core, 

(ii) mixing zone, and (iii) main region, as shown in Big. 2. 2. 
The potential core of the free jet is the region which 
extends from the nozzle exit to a position where the centre 
line velocit 3 ^ is about 93 percent of the velocity at the 
nozzle exit. Beyond this point the free jet centre line 
velocity decay is noticeable and the velocity decreases 
as the inverse of the distance. The theoretical velocity 
distribution of a turbulent round jet has been worked out 
by Tollmien [See Reference (104)] on the basis of Prandtl ' s 
mixing length theory [105]. G-ortler has also obtained the 
velocit 3 r profile of jet using Prandtl 1 s theorjr of free 
turbulence [See Reference (104)]. Ro analytical expression 
appears to be available for jet velocity in which the mixing 
layer and the main region are smoothly joined to each other. 
The literature survey reveals that most of the work is 
available on single-phase jet, but not much work has been 
done on the gas-solid particle jet. 

In the present work our interest is to obtain the 
centre line velocity of the solid particles in gas-solid 
particle jet. If centre lino velocity of the carrier phase 
(air) is known, then by using particle-phase momentum equa- 
tion, the centre line velocity of the solid particles can 
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be estimated. For this purpose an empirical equation 
proposed for evaluating fluid velocity at the centre line 
of a single -phase air jet has been used. This velocity 
is then modified to take into account the effect of admixture 
particles. In view of these a few additional assumptions 
have been made along with those indicated in Sec. 2.2. 

The concentration of admixture in the potential 
core remains constant while in the mixing zone and in the 
main region of the jet it varies due to entrainment of the 
surrounding fluid. The fluid phase velocity decays in the 
downstream of the jet due to transfer of momentum to the 
solid particles of the jet and due to the entrainment effect. 
The decrease in the fluid phase velocity due to transfer of 
momentum has been assumed to occur only in the potential 
core where the solid particle concentration is high. In 
the main region the concentration of solid particles is 
low, the decrease in fluid phase velocity has, therefore, 
been assumed to be mainly due to the entrainment effect. 

The concentration and the velocity profiles are assumed to 
be universal in nature. 

Albertson et al. [106] have given an empirical expre- 
ssion for the centre line velocity of axially symmetric, 
round air jet. Albertson's equations in a slightly modified 
form are 
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1 for x < x. 


(2.85) 


for x > x , 


( 2 . 86 ) 


where x is the distance from the nozzle tip, u is the 

7 gm 

centre line velocity of air in the jet, x is a constant 

c 

equal to the length of the potential core (x has dimen- 

c 

sion of length) and u is the velocity of air at the 

nozzle exit. When applied to gas-solid jet, the centre 

line velocity of the gas does not remain constant in the 

potential core (x < x ) but decreases due to momentum 

transfer to the solid particles, i.e., (u )_ 4 (u )^ . 

gm — o gm u \. — x 

c 

To take this effect into account, equation (2.86) has been 

modified by introducing another variable u^ in place of 

u . Here u represents the centre line velocity of the 

go go 

gas-phase at the end of potential core in a gas-solid parti- 
cle (two-phase) jet, i.e., (u ) ^ = u gc * The velocity 

u can be obtained by solving the total momentum equation 
S c 

of gas-solid particle mixture jet. Thus, equation (2.86) 
can be rewritten for the case of gas-solid particle mixture 
jet as 


for 


(2.87) 
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2.6.1 Potential Core Length 


Due to entrainment of the surrounding fluid, mixture 
density in the mixing zone and the main region of the jet 
changes considerably. Therefore, rather than using the 
term mixture ratio it is convenient to use the term concen- 
tration or mixture density. The concentration fi is defined 
as the ratio of mass flow rate of solid particle to the 
mass flow rate of fluid phase at any arbitrary section of 
the jet. Thus, the concentration 0 q at the nozzle exit is 
given by 




( 2 . 88 ) 


where a is the mixture ratio as defined earlier. 

To obtain the length of potential core in gas-solid 
particle mixture jet an approach similar to that used in 
single phase jet flow may be considered. Prom the theory 
of turbulent mixing proposed by Prandtl [105] it has been 
established that the thic.l-m.ess b of the turbulent mixing 
zone of submerged jet increases in proportion to the distance 
x from the mixing origin. Tho slope of the inner edge of 
the boundary layer of submerged piano jot with variable 
density is [104] 



2 + A 2 , 


(2.89) 
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where b-j is the thickness of the mixing zone and yj_ is the 
ordinate of the inner edge of the boundary layer at the end 
of the potential core. 'The quantities A , and A2 are 
expressed as 


g dp 


( 2 . 90 ) 


m 2 , 

; — gap , 


where 


(1 ~ p 1 * 5 ) 2 . 


TT Y 

here p = 


is the concentration profile. The term 


o mo 

is the ratio of the densities of the mixture at any arbitrary 

point in the boundary layer and in the core of constant 

p 

concentration ( potential core). This term ' a ~- can aUso be 

mo 

expressed as 


(1 + eO 

(1 + (2f 0 . 


( 2 . 91 ) 


Since p = 


, therefore. 
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A 


and 


A 1 = 


A„ = 


For round jet 


a _+_o.5^iy 
~'TT'+J 0 T ’ 


(0.45 + 0.129 0 J 


(1 + 0n) 




(0.316 + 0.067 0 ) 

rr+1 -y 


(2.92) 


7l 


JLO 

2 " 


(2.93) 


Substituting for A Q , and A ^ from equation (2.92) and 
for yq from equation (2.93), equation (2.89) gives 


b- 


d (1+0) 
no r o' 


2(0.416 + 0.309 0 O ) ' 


(2.94) 


It has also been shown semi-empirically fl04, 106] that 


(1 + 0.5 0 O ) 

! 5 itto - 


(2.95) 


where is an experimental constant which varies between 

0.2 to 0.3- 'The smaller values of c r are characteristic for 

5 

co-flo wing while the higher values are for counter flowing 
jets. 
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Substitution of equation (2.94) in equation (2.95) 
and taking c,- = 0.27 as obtained by Albertson et al. 1*106], 
the potential core length of a gas -solid particle mixture 
submerged jet is 


x 

c 


1.85 d n0 (1 + 0 Q ) : 


(1+0.5 0 O ) (0.416+0.309 0 O ) 


( 2 - 96 ) 


Substitution of 0 Q 


— ~ — r in equation 
(1-cO 


(2.96) gives 


1.85 d. 


x 


no 


(1-0.5 a) (0.416-0.107 a) 


(2.97) 


2.6.2 Flow Governing Equations of Gas-Solid 
Particle Free Jet 


In the jet region of inviscid gas-phase, neglecting 
body force, pressure gradient and taking = 1 for u < 
as discussed earlier, the conservation of momentum equation 
(2.40) for one- -dimensional steady state flow gives 


du. 


P u 
P '0 


dx 




du r 
P g u g dx 


- = 0 


( 2 . 98 ) 


Consider an imaginary stream tube of infinitesimal 

cross-sectional area in the centre of potential core 

c 

throughout its length. Assuming the concentration of 
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admixture in this stream tube to "be constant, the mass 
flow rate' of gas and solid particles through this tube 
are 


and. 


M 

P 


p ^ u 
P P 


AA C 


S 


> u 

S S 


AA 

c 


(2.99) 


( 2 . 100 ) 


The concentration 0 in the centre of the potential 

core is 




% p P \ 

m P n 

M g g g 


( 2 . 101 ) 


Using equations (2.88) and (2.101), equation (2.98) for the 
centre of potential core becomes 


du 


gm 


dx 


a 


dr 


^pm 


(1-cO dx 


for x < x 


- c 


( 2 . 102 ) 


where u^ m is the centre line velocity of particles in the 
jet. 

For one-dimensional steady state inviscid flow of 
gas -solid particle mixture jet in the centre of potential 
core , equat ion ( 2 . 41 ) gives 



85 


du pm 

dx 



u 


gm 


u 


pm 


1 ) 


(2.103) 


where is the time constant for momentum tr an sfer in 

du pm 

jet region. Substituting for d ’ rom equation (2.103) 

ir.t o e quat ion (2.102) gives 



, *gm _ . 

(1"«) 1 V " 


(2.104) 


”1 0~b 

■The time constant for momentum transfer in the jet 

region can be obtained using equations (2.7) and (2.10) as 


3 


jet 


1 

4 


i-V 


2 

P 


d * P. 


P 


[28 


(H 3et)-°. 85 


+ 0.48] I 


.jet 
Re * 


(2.105) 


Here particle Reynold’s number in jet region can be 

obtained by equation (2.5) in which o is replaced by o 

g J - ^ go 

which represents the material density of air at nozzle exit 
conditions i.e, at atmospheric conditions. ‘Thus, 


3 e "fc _ ^g o ( u gm “ u pm) 
i 'Re ~ * ” h" 


(2.106) 
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Substitution of equation (2.97) into equation (2.87) 
Aives the following equation for the centre line velocity 
of gas phase in the main region of the jet ; 


u 


gm 


1.85 d. 


'no 


(1-0.5 a) (0.416-0.107 a) 


'] 


u 


gc 


(2.107) 


for :c > T. . 

2.6.3 Flow Governing Equations in Dimensionless Form 
for Gas- Solid Particle Free Jet 

Equations ( 2.103) and (2.104) are the flow governing 
equations for the potential core region of the jet i.e. for 
x < x . while equations (2.103) and (2.107) are the flow 
governing equations for the main region of the jet i.e. for 
x > z . Equation (2.97) gives the length of potential core. 
For simplicity these equations will also he represented in 
dimensionless form vising dimensionless variables as given in 
o quat ion (2.64). 

Equation (2.103) rewritten in dimensionless form 

"becomes 



vf 



9 


(2.108) 
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where x* is the dimensionless distance from the nozzle 
exit, j.s a dimensionless constant, and u* m are the 

centre line velocity in dimensionless form for gas and 
particle phases respectively in the gas-solid particle 
mixture free jet. The dimensionless constant B* can be 
expressed as 


B 


jet 


no 


u 


gl 


1 

4 


p g 


no 


j et 

TvT'J 


Jet o . 8 5 


d P 


u 


{ Ee [28(JS£; )~ u * o:? +0.48] 


P gl 


(2.109) 


Similarly, equation (2.103) rewritten in dimensionless form 
becomes 



uB* 

(T7) 


u 


gn 


1) for 


u 


nrn 


•* * 
x < x c » 


(2.110) 


where x* is the potential core length in dimensionless form. 

The potential core length in dimensionless form can be 
expressed as 


„ r * 

'"c 


■v* 

J-i. 

c 



1.85 _ 

(1-0. 5a) (0.416-0. 107 k) 


( 2 . 111 ) 


Equation (2.107) rewritten in dimensionless form becomes 

u * = r , , 2 -- 3 5 _ . n 

m - (1-0. 5k) (0.416-0. 107a) 

*X" 

for x > x 
c 


9 


( 2 . 112 ) 
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where u* is the dimensionless centre line velocity of 
g° 

the gas phase at the end of potential core i.e. at x' v =x*. 

The particle Reynold’s number for jet centre line 
in terms of dimensionless variables is 


R 


jet 

Re 


U gl , * 

- SIB 


g 



(2.113) 


2.7 Numerical Solution Technique 


Equations (2.68), ( 2 . 69 ), (2.71), (2.73), (2.74), 

(2.75) and (2.84) are flow governing equations in dimem- 
sionless fora which describe the gas-solid particle mixture 
flow through a convergent-straight nozzle. Equations (2.108), 
(2.110) and (2.112) are the dimensionless flow governing 
equations for the gas-solid particle mixture free jet. These 
systems of nonlinear differential and algebraic equations 
can be solved numerically to obtain the particle velocity 
in the nozzle as well as the jet region. 


The gas phase velocity u p .^ and density at the 

nozzle inlet can be obtained from eqxiations (2.66) and 
(2.67), respectively. At the nozzle inlet the particle 
phase velocity u^ has been taken to be c Q times the carrier- 
phase velocity u gl . This factor c s may be called particle velocity 
ratio and defined as the ratio of particle-phase velocity to 
the carrier-phase velocity at the nozzle inlet. 
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In abrasive jet machining the solid particles are 
fed ' r ron rest (zero velocity) into the mixing chamber where 
it is mixed with the carrier gas. The mixture of gas and 
solid particles are carried through a tube or pipe to the 
nozzle. In doing so the solid particles also get accelerated 
by the carrier gas and attain some velocity at the nozzle 
inlet. The magnitu.de of the velocity attained depends on 
the distance travelled by the solid particles in reaching 
the nozzle inlet and the velocity of the fluid in the convey- 
ing tube. Hence the velocity ratio c will be a function 
of the length of the conveying tube and the carrier fluid 
velocity/ in the conveying tube. 

The numerical solution for evaluating the particle 
velocitjr can be carried out in two stages. In the first 
stage, the velocity of particle phase and gas phase is 
computed in the nozzle region. With the information at the 
nozzle exit (from the first stage of solution), the second 
stage of the solution gives the velocity in the jet region. 

The variables ifi, u*, [if , Tf , pf pf. and at the 

s P & P & u 

nozzle inlet are Imown and the increments au* aT* , Au c , 

V p» *V 

and ap* can be calculated vising equations (2.75), (2.84), 
(2.74), (2.71) and (2.72), respectively for the first interval 
of ax* using Runge- Putt a formula for first order simultaneous 
differential equations. Por each incremental step, the values 

di* 

of 31' , HR, A* and — 57 are computed from equations (2.76) to 
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(2. 79). The increments in the succeeding intervals can 
he calculated exactly in the same manner except that the 
end of the previous interval is taken as the starting point 
for the next interval. In this manner, the particle phase 
and the gas phase velocities at the nozzle exit can be 
obtained . 

In the second stage of the solution, equations (2.108) 
and (2.110) can be computed using Runge-llutta formula for 
first order simultaneous differential equations to get the 
particle phase and gas phase velocities upto the end of 
potential core i.e,, x* = x*. From this the gas. phase velo-- 
city u* at x* = x* is obtained. Using the computed value 

G 

of u* , the gas phase velocity at sections x*> x* in the 

£>C C 

main region of the jet can be calculated from equations 
(2.112). Knouing the gas phase velocity in the main region 
of the jet, the particle phase velocity for x* > x* can be 
obtained by solving equation (2.108) applying Runge-hutt a 
formula for first order different ial equation. 



CHAPTER 3 


THEORETICAL ANALYSIS OP MATERIAL REMOVAL RATE 
3.1 Introduction. 

Material removal in Abrasive Jet Machining occurs 
due to erosion caused Toy the impact of high velocity stream 
of abrasive particles onto the -workpiece. The erosion in 
such situations is a complex phenomenon and is not fully 
understood. Adequate data for governing parameters are 
not available. Some of these governing parameters are 
interdependent and are often difficult to control. Intui- 
tively, one might expect a large number of variables to 
enter the general description. Material properties (e.g. 
toughness, hardness, stiffness etc.), prior state of target 
surface (e.g. flaw density, level of surface stresses etc.), 
indentation parameters (e.g. indenter geometry, velocity of 
impact or rate of indentation etc.) are some of the factors 
which need to be considered. The complexities due to many 
interdependent variables seem to defy simplification for 
an experimentalist seeking to carefully separate variables 
and for a theorist attempting to accurately model the 
erosion phenomenon. Inspite of these difficulties, in the 
present work, an attempt has been made to develop a theore- 
tical model of erosion of brittle materials. Using the 
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developed erosion model an estimate of material removal 
rate during Abrasive Jet Machining of glass has been made. 

3.2 G-eneral Feature of Elastic-Plastic Indentation 
Fracture 

Although a brittle material does not show plastic 

def ormation, indentation fracture studies of brittle 

materials [48-51? 62,63] have shown that plastic deformation 

does occur even in very brittle materials like ceramics and 

glasses. Oh et al. [49] have shown that rather than fully 

elastic crack formation the elastic-plastic indentation 

occurs in brittle materials when the diameter of ind enter 

is below a critical value d (for soda lime glass 400 pm > 

pc 

d > 18 urn). Thus, under solid particle impacts the mode 
of fracture and deformation depend not only on the velocity 
and properties of the impacting particles relative to those 
of the target but also on the size of the impacting particles. 

Indentation studies [48,51,62,65] have shown that 
different sets of crack system are formed on the loading 
as well as unloading part of the indentation cycle. The 
sequence of events during indentation with a small spherical 
ball ind enter for one complete indentation cycle are shown 
schematically in Fig. 3»1» The contact of the sphere with 
the surface is initially elastic (Fig. 3*1 a) and is accu-- 
rateljr described by the Hertzian stress field. On increasing 
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the load, the stresses increase in proportion to the mean 
indentation pressure. Eventually the maximum shear stress 
exceeds the yield stress of the material and permanent 
deformation occurs (dark zone in Fig. 3.1 b). 

During the loading cycle, surface ring cracks occasio- 
nally form around the area of contact. On continued loading, 
these ring cracks sometimes develop into cone cracks but 
generally a zone of plastically deformed material develops 
around the indenter and a median crack develops within the 
target beneath the area of contact (Fig. 3.1c). This median 
crack also grows on continued loading. 

On unloading, the median crack begins to close-up 
owing to the compressive forces of the elastically deformed 
material around the indentation site (Fig. 3«1 d). However, 
on the surface, the elastic stresses are reversed and at 
about half the maximum load a system of radial cracks forms 
around the indentation site (Fig. 3.1 e). Most of the 
radial cracks are initiated from the edge of the area of 
contact and a few are extension of the subsurface median 
cracks. The system of radial cracks continue to propo’gate 
on unloading. The presence of Hertzian cone crack, if any, 
does not appear to influence the shape or extent of these 
radial cracks. As unloading proceeds the walls of the median 
cracks move closer but fracture debris and residual tension 
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once more prevent total closure. Indeed, the residual 
stresses may actually drive any partially formed median 
orach to full development. Just prior to complete un- 
loading the residual stresses (subsurface tension component) 
become sufficient to initiate and propogate a completely 
new sj^st era of sideways spreading saucer-like lateral cracks 
(Fig. 3* If). These lateral cracks continue their spread when the 
indent er load is finally removed, and may actually intersect 
the free surface to produce a chip [107]. These cracks 
form the basis of material removal during impact erosion 
of brittle materials. 

3.3 Plastic Deformation Threshold 

Erosion depends very much on the impacting particle 

velocity. If a small spherical particle (d < d ) impinges 

p ~ pc 

onto a surface of brittle material above a certain velocity, 
elastic-plastic deformation occurs and radial and lateral 
cracks system are formed. Thus, there exists a velocity 
threshold to cause plastic deformation and formation of 
cracks system. 

As soon as the spherical particle touches the surface 
of a flat body, stress concentrations are setup as a result 
of elastic deformation. Hertzian stress analysis indicates 
that the maximum stress concentration occurs in the centre 
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of the contact area at a depth of about half the radius of 
the projection of the contact area on the surface. As long 
as this maximum stress does not reach the yield strength 
of the target material only elastic deformation occurs and 
the collision is said to he perfectly elastic. On increase 
ing the load when the maximum stress in the centre of the 
contact area reaches the yield strength of the target 
material, the target starts deforming plastically. The 
impact -velocity at which the plastic deformation takes 
pla.ce can be related to the impact parameters using elastic 
collision formula of Hertz [108]. 

If the colliding bodies are spheres with one of the 
spheres having radius and mass infinitely large compared 
with that of the other, the equation expresses the colli- 
sion of a sphere with a flat body. In such a situation, 

Hertz equation gives the maximum contact force (F ) 

e max 

during elastic collision- as ( Appendix-II ) 


1 


(F ) 

' e max 


= 0 


20 % 


P ) 

P 


d‘ 



u. 


6 

5 

l pim 


(3.1) 


where u is ib® 'velocity, is the material density 

and d is the diameter of impacting particle. Here E and v 

P 



97 


are mod-ulus of elasticity and Poisson's ratio respectively. 
Suffix t represents the target material while p represents 
the particle material* 

Hertz has further shown that the radius a Q of the 
projection of the contact area on the body surface is 


a 


e 



(!' ) 
e max 






(9.2) 


as 


The average surface pressure q_ av can be expressed 


1 


av 


^e^max 
x a* 


(9.3) 


Substituting equat ions ( 3.1 ) and (3*2) in equation (3.3) 
gives 


40 p 


P 


av 


{7l(- 


V, 


E, 


1- 


+ 


i),4 


■A/5 2/5 

J u pim 


(3.4) 


Davies [109] has shown that the shear stress under 
a spherical indent er has a maximum value of 0.47 q. av at 
approximately 0.5 a below the surface of contact. At the 
onset of plastic deformation this shear stress must reach 
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the value of shear stress a ^ in simple tension at yielding 
Thus , 


= °- 47< l a v • (3.5) 

Usi^g liises criterion, Davies has further shown that the 
plastic deformation sets in if a ^ 0.5 Y + , which gives 

* = 1.06 Y. , (3.6) 


where Y^ is the yield stress of the target material in 
simple tension. 

Substituting the value of q from equation (3.6) 

a v 

into equation (3*4) and setting u . equal to u gives 

jP [pc 


u 


UC 


=■ 4.97 


5 

.2 

L t 
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P 


0* 


[' 


v: 


Dj 


1- 

4* 


Pi 2 


JCJ 


(3.7) 


P 


Hero u. is the threshold velocity at which plastic defor- 
pc 

mat ion in the target material is initiated when it is 
impacted by a spherical particle. The yield stress Yj_ 
can be related to the Vicker's hardness number . In 
particular Marsh r 110] has suggested that for glass the 
ratio of Vicker's hardness number to yield stress is approxi 
matcly 1.5 rather than 3 as found conventionally for metals 

H 

[ill]. Substituting Y, = , for glass target, equation 

(3.7) gives 
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u 


DC 


= 1.8 H. 


5/2 


(~r 


■)* r . 




1- y. 


E, 


E. 


£] 2 


(3.8) 


P 


3.4 Quasi-static Approach to Dynamic Impact 

Radial and lateral cracks formed during dynamic 
impacts [112] exhibit some similarity to that observed 
during quasi-static indentation [48,50,51,62]. These 
cracks are attributed to the development of an elastic- 
plastic stress field similar to that of a spherical cavity 
under internal pressure [50] with cracks propogating on 
planes normal to the maximum principal tensile stress (in 
the same manner as the Hertzian cracks are formed during 
the elastic indentation). The lateral cracks usually 
initiate and develop during the unloading half cycle, and 
are motivated by the residual elastic-plastic stress field. 
Further, the lateral cracks are sub-surface cracks and run 
parallel to the free surface below the impact site. The 
lateral cracking is more difficult to analyse because it 
is subsurface and primarily dependent upon the proximity 
of the freo surface. Due to these complexities in the 
elastic-plastic residual stress analysis, the propagation 
of the lateral cracks has not yet been fully explained, 
[50,51]. Under impact conditions the stress-fracture inter 
relations in the plastic response regime are expected to 



100 


Toe even more complex than for the quasi-st at ic indentation. 
Further, the shock waves, plastic and elastic waves and 
their interactions with rapidly moving cracks during impact, 
complicates the analysis. Becatise of these complexities 
the quasi-static indentation approach has heen used to pre- 
dict the erosion loss under impact situations in the present 
work. This approach is similar to that used by Evans and 
Wilshaw [50]. It is assumed that lateral cracks gjrow in a 
quasi-static manner as a result of residual stresses produced 
by the impact event. The elastic Herts equations are 
extended for plastic load also. This approximation has been 
found to be reasonably valid for prediction of erosion loss 
at low velocity ( <160 m/sec.) impacts [112], 

3.5 Mechanics of Erosion 

During solid particle impact erosion, the material 
removal may occur due to crushing of target surface Linder 
impact and/or due to interaction of lateral cracks. Crush- 
ing part of erosion is more dominant for large particle 
impacts, while the lateral crack interaction is mainly 
responsible for erosion loss for small particle (d^ < <3-p C ) 
impacts. Under such situations when impacting particles 

are small (d < d ) the erosion volume can be obtained by 
P ~ P c 

knowing the depth and size of the lateral crack. 
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3.5.1 lateral Crack Depth 

Consider an elastic-plastic collision where the 
colliding particle is spherical in shape and deforms only 
elastically, while the target surface deforms both elasti- 
cally and plastically. For elastic-plastic impact, using 
Hertz analysis Bitter [16] has shown that the depth of 
elastic penetration z 3 (Fig. 3.2) is related to the maximum 
surface pressure occurring in the centre of the contact 

area and is given by 


z„ = 


TC 

8 


( 1 - 


»t> 


(1- 


L max 


Ej 


E 


i > ]2 


(3.9) 


According to Hertz the average surface pressure 
2 


l a v — 3 ^max" 


The largest specific load that can be 
exerted upon a body without catising plastic deformation is 


called the elastic load limit ( a , ) . Therefore the highest 

yt 

average surface pressure that can be reached during a 

p 

perfectly elastic impact is -- a , . Thus, the elastic load 

3 yc 

limit can be expressed in terms of average surface pressure 
as 



av 


1.5 q., 


( 3 . 10 ) 
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Substituting for ci from equation (3.6), equation (3.10) 

av 

gives 

°yt = 1-59 T t . (5.11) 

It is assumed that the target material undergoes 
plastic deformation at constant yield stress without strain 
hardening. After elastic load limit is reached, plastic 
deformation sets in at the centre of contact area. On 
further penetration of the particle into the workpiece, 
the area which is loaded to the stress value c, increases. 

y ° 

The total depth during elastic -plastic impact will, 
therefore, be 


z t = 


z Q + z - n 
e p 


( 3 . 12 ) 


where z^ is the depth of concavity due to plastic deforma- 


z ion. 


At the end of loading cycle during impact, the edge 
of the concavity (edge of the plastically deformed area) 
is just loaded to a pressure equal to cr Outside the 
plastically deformed contact zone the sphere also touches 
the target surface, but the specific load decreases from 
Oy^_ at the edge of the concavity to zero at the edge of 
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the total contact area and the deformation is purely 
elastic (fig. 3*2). 

Let the projection of the total contact area normal 
to the target surface has a radius . In practice z and 

G 0 

z are small compared to the particle size d . Therefore, 
p P 

the total contact area may be approximated by 

K md p (z e + z p ) . (3.13) 

Similarly if the radius of the projection of the plasti- 
cally deformed area is a , then 

p 

ma p jw xd p z p . (3.14) 

The elastically loaded area can be obtained as 


Ia t 


7i a 


V 


% d_„ 

P ^ 


(3.15) 


Equation (3-15) shows that during elastic-plastic 
impacts the size of the contact area subjected to elastic 
deformation only is independent of the size of the inden- 
tation formed. Further, this is equal to the Hertzian 
contact area in a perfectly elastic impact in which 
just equals If we assume the average indentation 
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2 

pressure q_ av in this area to be equal again to ^ cr^, 
equation (3.4) can be used to calculate the elastic impact 
velocity u and thereby the elastic part of the energy X g 
absorbed in this region. The elastic impact velocity can 
now be expressed as 


u 


pe 


TC 

2fl0 


3Tt 



( 1 - 




( 1 - 


B, 


4 * 


E. 


.D] 2 


P 


(3.16) 


Substitution of equation (3.11) into equation (3.16) 
gives 




(3.17) 


Thus, the impact energy spent in elastic deformation 2 is 

X e = i \ > (5- 18 ) 


where mass of impacting particle = -g p^d 


P P 


The elastic reaction forces are set up at the end 
of the impact cycle as a result of this elastic impact 
energy. The elastic reaction forces then press the particle 
out of the body (target) giving the -"particle an amount 
of kinetic energy equal to the stored elastic deformation 
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energy. From the energy balance of collision cycle, 

Bitter Tl6] has been able to show that the energy absorbed 
by the target during plastic deformation X^ is 


p ~ 


( n . 

• im 


fx e ) 


(3.19) 


where X ±m is the impact energy of the colliding particle. 

Formation of plastic indentation requires an amount 
of energy equal to the plastic deformation work in penet rat- 
ing the particle to a depth ( z p) maz during indentation cycle. 
This gives 



^ z p^max 

/ % V az 


7td ( z ) 2 a , 
p p 'max yt 


( 3 . 20 ) 


If the particle impacts onto the target surface with 
a velocity u^, the impact energy of the colliding particle 
is 


X. 


im 


i 

T 



( 3 . 21 ) 


Substitution of X^ from equation (3.18) , X from 
equation (3.20) and X^ from equiation (3. 21) into equation 
(3.19) gives 
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0 

^p^max = ^ dp ^ Upim ~ 


( 3 . 22 ) 


Using equation (3.14), the maximum radius of projection 
of the -plastically deformed area, (&_)„„_ can he expressed 

p JUSL.A. 

as 


• a p^max 


^p^ z p^max 


(3.23) 


Substituting for (z_.)„,,_ from equation (3.22) into equation 

p max 

(3.23) gives 
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(3.24) 


Hill's analysis [113] for the expansion of a spherical 
cavity under internal pressure in an elastic-plastic solid 
was extended by Johnson [114] to analyse the elastic -plastic 
indentation. In Johnson’s analysis the analogy of expansion 
of a spherical cavity with internal pressure has been taken 
for the plastically deformed zone (hydrostatic zone) about 
the indent er. According to Hill [113] the size of the 
plastic zone C Q is expressed as 


E, 


Si. 
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-3(1- ^lxT ] 
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(3.25) 
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Outside the plastic zone the material is loaded only elasti- 
cally. Thus, the elastic-plastic boundary occurs at a depth 

C below the target surface, 
o 

Ev an s et al. [63] and others [50,51,53] have suggested 
that the lateral cracks initiate at the elastic-plastic 
boundary. Thus, the lateral crack depth can be approxi- 
mated bir the size of elastic-plastic boundary which gives 


E. 
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5Tl- 
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( 3 . 26 ) 


Prom the available data of mechanical properties for 

n 

,qi/3 


oiTy."J 


i 

different glasses, Table 3.1, the quantity [-yz — - 

J V 1"' ^4- J 

has been calculated and found to be constant ( ^ 2). Thus 

the lateral crack depth for glasses and similar brittle 
materials (ceramics) can be expressed as 


Z 1 ~ 


2 a.. 


(3.27) 


For full penetration at the end of impact cycle a = (a ) . 

_P ITlSlX! 

Substituting for (a )_, „ from equation (3.24) into equation 

j) in six 

(3.27) gives 


^> 1/4 ( 


y* 


, u . -~u ) : 

p pun pc' 


(3.28) 
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Substituting for a , from equation (3.11) into equation 

7 H t 

(3.28) and taking the expression for lateral 

crack depth into the glass target is 


z 
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1.26 t /)^ 4 d . 


(u. 


-u 


pom pc 


)■ 


(3.29) 


Equation (3.29) gives the estimate of lateral crack 
depth in glasses and similar brittle materials (ceramics). 
This expression is similar to that obtained by Evans et al. 
[63] » with the only difference that the present expression 
takes into account the velocity threshold effect. Equation 
(3-29) has been found to be in good agreement with the 
experimental data for lateral crack depth (Table 3.2). 

3.5.2 lateral Crack Length 

During unloading part of the indentation cycle 
irreversibly deformed material within the plastic zone is 
directly responsible for the formation of a residrial elastic 
stress field which provides the driving force for lateral 
crack extension. Any numerical or rigorous analytical 
solution of the stress fields due to plastic indentation 
is likely to be complicated by non-linear deformation 
behaviour. The complexity of anals^sis is further increased 
because of the dynamic nature of the problem and the effect 
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of indent er geometry on the material behaviour. In view 
of these complexities and the current state of knowledge 
of the indentation stress field, a semi- empirical approach 
has been adopted to express the lateral crack length in 
terms of impact parameters, 

Evans and Wilshaw [50] have expressed the normalized 
lateral crack length in terms of dimensionless functions 
of different variables influencing the crack extension 
phenomenon as 


f 


C-, 

(~i) 

1 a p> 




f "' 1 ( 3 > t ) 


(5.30) 


where f are dimensionless functions, X is stress inten 
sity factor,^ is coefficient of friction at the indenter 
and target interface and y, is Poisson's ratio of the 

"G 

target material. 


The uniaxial yield stress is related to the hardness 
of the material by [115] 


H t 3 3(E,/Y,) 

= 0.28+0.6 [ — - — — — — ] ln[ — — - — — 

t 3»(l-2 » t )(g^) (4+ a> t )-(^)(l-2^)(l +v ± ) 


(3.31) 
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- 1 * 

which for small |p ( < O.l) reduces to 



E, 


0.6 ln(~) +0.28 + 0.6 ln(^^-) 
t t 


Q 


(3.32) 


Here q is a material constant. Substituting for Y, from 

u 

equation (3.32) and setting K equal to K (the crack 
propogation condition) equation (3*30) gives 


1 P 


K Q 
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H t^ 


It, 
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(~°-) f, 

a P 5 


■i 


b t ) 


■4 


CO]. 


(3.33) 


Here Y, q is the fracture toughness of the target material 
which represents the resistance to crack propogation and 
is a material property. In this equation the terms within 
the square bracket remain almost constant since the friction 
coefficient r for a fixed profile ind enter is only marginally 
influenced by the material, properties and may be assumed to 
be constant, Poisson’s ratio being small (0.22 - 0.3) 

does not exert a dominant influence on the residual stress 

C 

field, and the term is essentially a constant. Evans 

a 'p 

and Wilshaw ^ 50 ] have ta,kon functions f£, and ^4 ©Ra- 
tion (3.33) as approximat ely constants and independent of 
the material and have expressed tho normalized lateral crack 
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length as a function depending primarily on the quantity 
K 0 

—S ~ — They have obtained a plot of the normalized lateral 
^t ^P C n 

crack length — • as a function of the normalized fracture 
a P 

toughness If Q/H.YaL for a variety of target materials as 
shown in Pig. 3.3(h). They have further shown that the 
lateral fracture is independent of the crack formation load 
and the surface condition (surface flaws) hut dependent on 
the plastic zone and the residual stress in the vicinity 
of the elastic-plastic boundary. Once plastic indentation 
has occurred the peak stress and the stress distribution 
remains essentially constant (except for the increase due 
to work hardening), and only spatial extension of the 
stress field takes place as indentation proceeds. Thus, 
the applied load (contact force) can influence ^q(^p/ a - 0 ) 
value only through the spatial stress field extension. 

This specifies that the probability of crack propogation 
increases as the volume of stressed material as determined 
by the impression diameter, increases or when the critical 
crack propogation stress decreases. The volume dependence 
thus provides a critical force condition, and the critical 
stress requirement provides a micro structure dependence. 
Using the plot of Pigs. 3.3 (a & b) and relating the target 
hardness H,_ with the contact force P during plastic inden- 


tation, for 


2, Evans and Wilshaw [50] have expressed 
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the lateral crack length by the following approximate 
proportionality : 

0 ± cc ( i ) 5/4 . (5.34) 

Equation (3* 34) can be used to obtain pertinent information 
about the lateral crack extension if the contact force 
during plastic indentation is known. 

The contact force during plastic indentation does 

not remain constant but varies with the depth of penetration. 

Ho analytical solution appears to be available which relates 

E with z . Using numerical solution and fitting an empiri- 
P P 

cal relation for the contact force (non-work hardening target 
material), Evans and Wilshaw [50] have expressed the relation- 
ship between E and z as 
P P 

E =. 2,545 G, \T z 3 J 2 , (3.35) 

P c P P 


where G-^ is the shear modulus of the target. 

Eor full penetration during plastic indentation 
i.e., z = (z ) t equations (3.22) and (3*35) give, the 

JJ ILlct^X 

maximum contact force (F ) as 

p max 
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(F ) 

v p'max 


0.664 G + if. (— ) 

c ~ D v 




(u 


. —u ) 
pun pc 


3/2 


(3.36) 


Substituting F = (F ) , equations (3.36) , (5.34) 

JJ JJ Ulci-^. 

and (3.11) give 


* (?4) 3/4 Y /2 (|£)9/ 16 ( 

P X 4- 


U . -11 ) 

pun pc 


9/8 


(3.37) 


Substituting for in terms of the target hardness (Y^ 
for glasses), equation (3.37) can be rewritten as 
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<X ( K ) 


2 i p i_ i 

( U ) 16 (U . -u ) 8 


i ' 4 a 2 


P H, 


pim pc 


(3.38) 


3.6 Material Removal Rate 

Material removal during solid particle impact may 
be assumed to occur due to intersection of lateral cracks 
from adjacent impacts. The material is removed in the form 
of micro-size chips. These micro-size chips may be assumed 
to be cylindrical in shape having radius equal to the lateral 
crack length and depth equal to the lateral crack depth z 1 
produced under the impact- Thus, the volume of material 
removed per impact (v) may be expressed as 

= « C* z ± . 


V 


(3-39) 
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Substituting for z 1 and C 1 from equations (3.29) and (3*38) 
into equation (3.39) gives 
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3, 
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cC 


( if' V C ( 


11 

) 8 ( 


U. 


nun 


-u ) 
pc 


11 

4 


(3.40) 


Equation (3.40) gives the volume of material removed from 
the target siirface during a single impact. 

As mentioned earlier, in AJK high velocity air- 
abrasive mixture jet issued from a nozzle impinges onto 
the workpiece surface. Material removal occurs du.e to 
erosion caused by multiple particle impacts. Assuming that 
all particles in the jet during impact have velocity u 

pm 

equal to that at the centre of jet and have negligible 
particle-particle interaction effect, the cumulative effect 
of multiple particle impact can be obtained by direct multi 
plication of the number of impacts with the volume removed 
by a single impact . 

The number of impacts per unit time can be estimated 
from the number of particles passing per unit time at any 
section of the air-abrasive mixture jet. Thus, the number 
of impacts per unit time depends on the mass flow rate of 
abrasive particles in the jet, material density of the 
particle and size of the particle. The number of impacts 
per unit time can be expressed as 
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6 M P 

P n ~P 
P P 


(5.41) 


Considering multiple particle impact, the volumetric 
material removal rate in ATM can he expressed as 


V = 11 v 

r p 


(5.42) 


Taking u . = u , equations (5.40), (5.41) and 

pim pm 


(5.42) give 


V ot - 
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k d G 2 P 8 (u ”U ) 
p p t p pm pc 


11 




11 / 8 " 


(5.45) 


Equation (5.45) gives material removal rate for 
machining of brittle materials such as glass, ceramics etc. 
in terms of proportionality equation. When applied to real 
situations of abrasive jet machining process the effective 
number of impacts onto the workpiece surface per unit time 
is likely to be less than the number lST calculated from 
equation (5.41) because of the particle-particle interaction 
effect. The particle-particle interaction effect primarily 
depends on the stand-off distance. At small stand-off dist- 
ances due to less flaring and high particle concentration in 
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the jet the rebounded particles from the workpiece surface 
collide with the incoming particles of the jet and reduces 
the effective number of impacts onto the workpiece surface 
per unit time. Thus, the effective number of impacts will 
be negligibly small at nearly zero stand-off distance. When 
the stand-off distance is high the particle-particle inter- 
action effect is low gives large number of effective impacts 
approaching to the value E . Further, the particle-particle 

p 

interaction effect will be more at high mixture ratio due to 
large number of particles per unit volume of the mixture. 

But for the sake of simplicity the effect of mixture ratio 
on particle-particle interaction effect has not been consi- 
dered in the present analysis. Introducing a term X as a 
function of stand-off distance to take into account these 
effects in real situations, equation (3*43) can be rewritten 
as 
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2 2 

X M d G- 2 P 8 (u 
P P t p ^ pm 
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ii 
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( 3 . 44 ) 


The term X can be evaluated from a typical experimental 
data. 
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Pig. 3.2 Pressure distribution at contact area 

(«) elastic impact (b) elastic-plastic impact. 
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CHAPTER 4 


EXPERIMENTS 


4 . 1 Int r oduct ion 

Experiments have been conducted on an experimental 
abrasive jet machine specially designed and fabricated for 
the purpose. The present chapter deals with the design 
aspects of this abrasive jet machine and the experimental 
investigations of various input parameters affecting erosion 
"rate and nozzle wear. During the first phase of the experi- 
mental study, the effect of various input parameters such 
as stand-off distance, grain size, mixture ratio, carrier 
fluid pressure etc. on volumetric erosion rate, penetration 
rats, diameter of eroded cavity etc. have been studied. In 
the second phase, experiments have been carried out to 
investigate the effect of nozzle diameter, nozzle length 
and entrance angle on nozzle wear. The other parameters 
which have been considered for the wear study are mixture 
ratio and abrasive particle size. 

4.2 Experimental Set-up 

The schematic layout of the abrasive jet machining 
set-up is shown in Pig. 4.1. Dry air from the compressor 
(l) through dehumidifier and filter (4) is allowed to enter 
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the mixing chamber (8) through a pressure regulator (3). 

The pressure regulator enables regulation of inlet pressure 
to the system* The flow meter (6) and. the pressure gauge 
mounted on the outlet end of the pressure regulator indicate 
the flow rate and the pressure of air entering into the 
mixing chamber. Fine grained abrasive particles stacked in 
the main abrasive chamber (18) flow down under gravity 
action into the secondary abrasive chamber (17) through a 
metering tube (16). The metering tube maintains a constant 
head of abrasive particles in the secondary abrasive chamber 
irrespective of the abrasive particle head (stock of abrasive 
particles) in the main abrasive chamber. Abrasive particles 
from the secondary abrasive chamber flow through the abrasive 
metering venturi (15) > under the action of a mechanical 
vibrator (7), into the mixing chamber where it is mixed with 
the incoming air. The abrasive flow from the secondary 
abrasive chamber is controlled by adjusting the amplitude 
and frequency of the mechanical vibrator attached to the 
secondary abrasive chamber. The amplitude of vibration can 
be adjusted by adjusting the throw of eccentric driving the 
slider block of the vibrator unit. The frequency can be 
varied by varying the rpm of the motor driving the eccentric 
unit of the mechanical vibrator. The mixing chamber, the 
secondary abrasive chamber and the main abrasive chamber are 
maintained at uniform pressure through pressure equalizer 
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tubes. Hence, the abrasive flow in the system is not 
influenced by the pressure and flow rate of air and a 
precisely controlled air-abrasive mixture can be obtained 
at any desired nozzle pressure. The air-abrasive mixture 
is then directed through a flexible tube to the nozzle 
fitted in the nozzle holder. The nozzle holder facilitates 
easy, accurate and leak proof holding of nozzles. The 
mixture of air-abrasive particles appears in the form of a 
stream of air-abrasive jet at the nozzle exit. The nozzle 
and its holder (14) is mounted on vertical and horizontal 
slides for accurate location of the nozzle tip from the 
workpiece. The pressure gauge (9) connected to the nozzle 
holder reads the mixture pressure at nozzle inlet. The 
workpiece (ll) is mounted on a fixture (12) inside a closed 
working chamber (10). The working chamber is connected to 
a vacuum cleaner (13) for removal of abrasive along with the 
swarf. Photographic views of the set-up are shown in Pig. 

4 . 2 ( a & b). 


4.3 Design Considerations 

In AJM , suitably proportioned air- abrasive mixture 
from the mixing chamber is carried pneumatically to the 
nozzle through a flexible tube. For pneumatic conveying of 
solid particles through a pipe a certain minimum velocity 
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of the carrier fluid is required to keep the particles in 
suspension [115-118] • This minimum velocity is also known 
as the minimum transport velocity or the saltation velocity. 
When velocity of carrier fluid in any part of the transport 
system falls below the saltation velocity, the solid parti- 
cles nay settle out of the flowing stream and may cause 
irregular flow of mixture. Thus, the units of AIM through 
which air™abrasive mixture is flowing must be designed to 
maintain the saltation velocity for the specified values of 
mixture ratio, working pressure, type of abrasive particles, 
size of abrasive particles and the exit diameter of the 
nozzle. 

The saltation velocity for pneumatic conveying of 
solid particles in horizontal tubes can be estimated by 
the method proposed by Zenz [115]. In AJM under real situa- 
tions flow of mixture in conveying tube neither conforms to 
horizontal conveying nor to vertical conveying, but to a 
close approximation it may be treated as conveying through 
inclined tubes. According to Zenz and Othmer [118] there 
is no marked change in the saltation velocity when the 
conveying tube is inclined as much as 10 degrees from hori- 
zontal. Further, they have shown that for intermediate 
positions 22 to 45 degrees from the horizontal, the saltation 
velocity may be taken to be 1.5 to 3 m/sec higher than in 
the horizontal conveying. 
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In AJM the mixture ratios used are generally below 

0.5, the abrasive particles used are usually in the range 

of 10 to 50 jam and the carrier fluid pressure used is gene 

rally in the range of 196.2 to 588.6 MT/m 2 (gauge) [5-6]. 

In the present case the set-up has been designed for system 

inlet pressure of 0.5 Mt'T/m 2 (abs) with the nozzle exit at 

atmospheric pressa're [0.101 MF/m 2 (abs)]. For adiabatic 

flow through convergent nozzles it has been shown that when 

the nozzle exit pressure reaches a value equal to or less 

than 0.528 times the inlet pressure, the nozzle exit velocity 

becomes sonic (540 m/sec) [119 ]. In the present set-up 

convergent -straight nozzles are used with exit pressure of 

0.202 times the inlet pressure and hence, the carrier fluid 

will attain sonic velocity at the nozzle exit. Taking the 

mixture ratio to be 0.4, the maximum size of Al^O^ particles 

to be 50 jam, the minimum nozzle diameter to be 0.7 mm and 

the carrier fluid velocity at nozzle exit to be 540 m/sec, 

the saltation velocity u for horizontal conveying is 0.62 

gs 

m/soc as shown in Appendix-Ill. The minimum safe velocity 
of carrier fluid u..^ in the mixture conveying tube can, 
therefore, be taken as 5.62 m/sec. 

As mentioned earlier, the abrasive particles are 
fed from the secondary abrasive chamber, through the metering 
venturi which is connected to it, into the mixing chamber 
under the action of the mechanical vibrator. The mixing 



127 


c ha mber is held in the vertical position and is connected 
to the secondary abrasive chamber through the metering 
venturi. High pressure air enters at the top of the mixing 
chamber and is mixed with abrasive particles entering 
through the metering venturi which is connected on one side 
of the mixing chamber. The vertical orientation of mixing 
chamber facilitates in easy fall of abrasive particles 
under gravity action in the mixing chamber from the metering 
venturi. Thus, the flow of air-abrasive mixture in the 
mixing chamber is in vertically downward direction. 


The minimum velocity of carrier fluid for vertical 
conveying, known as choking velocity (u-g C - n ) , is generally 
taken to be equal to u o . g for conveying of uniform size 
particles, and 0.2 to 0.33 times u for conveying of mixed 

gs 

size particles [120]. Thus, for the working conditions 
stated earlier, the minimum safe velocity of carrier fluid 
n^ch in the mixing chamber has been taken to be 0.62 m/sec. 
On the basis of minimum safe velocities of carrier fluid in 
the mixing chamber and in the mixture conveying tube for a 
minimum nozzle diameter of 0.7 mm, the maximum diameters of 
the mixing chamber and the mixture conveying tube comeout 
to bo 7.3 mm and 3 mm, respectively as shown in Appendix- IV. 


The air compressor capacity is selected on the basis 
of flow requirement for the maximum nozzle diameter in the 
working pressure range. Since working pressure range for 
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carrying out the experiments has been chosen to be 344.863 
to 372.455 MT/m 2 (gauge) and maximum nozzle diameter as 
1.62 mm, a compressor with working pressure of 400 kET/m 2 
(gauge) and discharge capacity of 0.725x10" 3 m 3 /sec will 
be sufficient for the purpose. 

4.4 Preliminary Experiments 

The set-up was designed and fabricated taking into 
account the design considerations indicated and was initially 
operated for a few hours to check for leakage at various 
joints. After ensuring no leakage condition, the sot-up 
was used for experiments. 

In the present set-up, various values of abrasive 
flow rates were obtained by varying the frequency of vibra- 
tion. The necessary calibration to this effect was carried 
out. Por obtaining the flow rate, the abrasive jet was 
directed into water for two minutes at the specified fre- 
quency of vibration and then the abrasive particles were 
filtered and dried to obtain the weight. Calibration curves 

for abrasive flow rate for different particle sizes are shown 

/ 

in Pig. 4.,0. 

Preliminary experiments have boon carried out using 
^l 2 0 3 abrasive particles of grain sizes ranging from 25 to 
48 (am for both volumetric erosion as well as nozzle wear 
studios. The most suitable range of grain size appeared to 
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"be 30 to 38 [am. Coarse grains caused excessive damage to 
the mixture conveying tube while fine particles often caused 
choking. Though a wide range of mixture ratio (0.052 to 
0.508) has been used but most of the experiments have been 
carried out only in the range of 0.052 to 0.268 to avoid 
excessive damage to - the mixture conveying tube and for 
convenience in setting the motor rpm and the air flow rate. 
The procedure for calculation of mixture ratio is given in 
Append ix-V. 

All machining tests have been performed on ordinary 
glass tost pieces by directing the air-abrasive jet at 90 
degreees to the work surface using tungsten carbide nozzles 
of circular cross-section having 0.712 mm exit diameter and 
convergent -straight geometry (Pig. I.l). The machining 
(cavity drilling) tests have been carried out at different 
stand-off distances ranging from 2 to 20 mm. Beyond 20 mm 
stand-off distance, the top diameter of eroded cavity was 
not clear due to excessive flaring of the abrasive jet. 

In view of this all tests have been carried out upto a 
stand-off distance of 20 mm. For ease in measuring the 
cavity diameter and depth, each test was carried out for 
60 seconds. The volume of material eroded, the cavity top 
diameter and eroded depth at the centre have been evaluated 
by projecting the eroded cavity on the screen of a measuring 
pro j 0 ct or . 
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Nozzle wear tests have "been performed on high carbon 
high chromium tool steel nozzles of convergent -straight geo- 
metry (Fig. 4.4a). Experiments have been carried out to 
ascertain the wear rate of nozzle as well as wear profile. 
For this purpose a set of split nozzles were used. These 
nozzles were made in two halves and assembled together to 
get the final profile. The nozzles were disassembled after 
use and the inside surface was examined under a profile 
projector to obtain the wear profile. During preliminary 
experiments, the nozzle profiles were traced after three 
minutes intervals. Nozzle wear profiles obtained after 
various operating times are shown in Fig. 4.4 (b-e). It 
is clear that the profile of the worn out nozzle closely 
resembles a truncated cone. The wear volume in all subse- 
quent experiments was, therefore, calculated on the basis 
of the inlet and exit diameters of the flow channel by 
assuming a straight flare from the entry to the exit side. 

The wear volume was also obtained by weighing the nozzle 
before and after the use. The wear volumes thus obtained 
are plotted in Fig. 4.5. The wear volumes evaluated on the 
basis of the inlet and exit diameters of the nozzle flow 
channel are also plotted in this figure. It is clear that 
both the results closely match. 
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Nozzle wear profiles shown in Fig. 4.4(b~e) clearly 
indicate that the exit diameter of the nozzle flow channel 
becomes almost 1.3 times the initial exit diameter in approxi 
mately 20 minutes. Also the assumption that the wear profile 
resembles a truncated cone does not appear to be valid after 
15 minutes of use., All subsequent experiments were, there- 
fore, carried out upto 15 minutes only. Further, the time 
interval for recording of successive wear profile was also 
reduced to 2.5 minutes. 

4.5 Experimental Conditions 

After obtaining satisfactory results in the prelimi- 
nary experiments, erosion rate and nozzle wear experiments 
have been carried out under the following experimental 
conditions : 

ITozzle 

Material s Tungsten carbide and high carbon high 
chromium tool steel ( 2.15 % C, 12$ Or, 

0.35$ Mn, 0.2$ Si, 0.8$ Mo); 


Diameter 


0.712, 0.78, 0.8, 0.9, 0.96, 0.97, 1.17 
and 1.62 mm . 
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Carrier Fluid 

Carrier gas : Dry air ; 

Ini art pressure: 344.863 to 372.453 kh/m 2 (gauge); 
Nozzle pressure: 98.1 to 234*3 kN/m 2 (gauge). 

Abrasive 

Type i Aluminium oxide; 

Size : 25, 30, 38 and 48 pm; * 

Mixture ratio ; 0.C52 to 0.308. 

Yibrat or char act erist ics 

Amplitude s 1.9 mm; 

Frequency : 5 to 8 c/sec. 

Work Material : Glass. 

4*6 Experimental Results 

Experiments have been conducted to investigate the 
effects of various input parameters on volumetric erosion 
rate, penetration rate, top diameter of eroded cavity etc. 
and on wear of nozzle of various geometries. Experiments 
have been carried out over the range of experimental condi- 
tions stated earlier. All erosion tests have been carried 
out at system inlet pressure of 344.863 klT/m 2 (gauge) while 
nozzle wear tests have been conducted at system inlet pressure 
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of 372.453 Hl/rn 2 (gauge). The effect of stand-off distance, 
mixture ratio, abrasive particle size, nozzle pressure etc. 
on volumetric material removal rate, penetration rate, top 
diameter of eroded cavit^r are shown in Figures 4.6 to 4.14 
while wear results obtained with nozzles of various sizes 
are plotted in Figures 4.15 to 4.23* 

4.6.1 Erosion 

The effects of stand-off distance on volumetric 
material removal rate and penetration rate are shown in 
Fig. 4.6. The profiles of the machined cavity at several 
stand-off distances are also indicated. Figure shows that 
both material removal rate and penetration rate first 
increase with increase in stand-off distance and then 
decrease giving an optimum value. The maxima of the pene- 
tration rate and the volumetric material removal rate do 
not occur at the same stand-off distance for the specified 
machining condition. Figure further shows that the maximum 
penetration rate occurs at smaller stand-off distance than 
the value where maximum material removal rate is obtained. 
Photographs of the machined cavity at different stand-off 
distances are shown in Fig. 4.7. . 

Volumetric material removal rate not only depends on 
stand-off distance but also depends on the mixture ratio. 

Fig. 4.8 shows the variation of volumetric material removal 
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rate with, stand-off distance for different mixture ratios. 

I'c shows that the optimum stand-off distance at which 
material, removal rate is maximum does not change much with 
variation in mixture ratio. 

The variation in penetration rate with stand-off 
distance for different mixture ratios (fig. 4.9) indicates 
that the penetration rate increases with increase in mixture 
ratio. The stand-off distance corresponding to the maximum 
penetration rate also increases with increase in mixture 
ratio . 

The top diameter of the eroded cavity also varies 
with stand-off distance and mixture ratio as indicated in 
Fig. 4.10. It is clear that the top diameter of the eroded 
cavity increases with increase in mixture ratio and stand- 
off distance but tends to saturate beyond a certain stand- 
off distance. For low mixture ratio the saturation tends 
to occur at smaller stand-off distance than for high mixture 
ratio . 

The erosipn rate is very much dependent on mixture 
ratio and particle size. To have a comparative effect of 
abrasive particle size on volumetric erosion rate, the volume 
erosion factor e was employed. This is defined as, 
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V olume tric material removal rate (mm 3 /min ) 

e v ~ Weight of abrasive particle impacted per -unit time(g/min) 

(4.1) 

Pig. 4.11 shows the variation in this erosion factor e with 
stand-off distance for different grain sizes. It shows that 
the volume erosion factor increases with increase in stand- 
off distance and then falls with further increase giving an 
optimum value. Further, this optimum stand-off distance 
does not appear to vary when grain size is varied. The 
coarse grains, however, give higher values ox erosion factor 

£ . 

V 

The variation in optimum material removal rate with 
mixture ratio for different grain sizes (Fig. 4.12) indicates 
a linear increase in material removal rate upto a certain 
value of mixture ratio and then tends to fall. Again coarser 
grains give higher value of volumetric material removal rate 
while finer particles have larger linearity range. 

Fig. 4.13 indicates the effects of grain size on 
penetration rate. It is clear that coarse particles give 
higher penetration rate. The stand-off distance where maxi- 
mum penetration rate is achieved, however, shows only margi- 
nal effect of grain size. 



136 


The effect of nozzle pressure on volumetric material 
removal rate is shown in Pig. 4.14. The rate of increase 
of material removal rate appears to decrease with increasing 
nozzle pressure and has a tendency to saturate beyond a 
certain value of nozzle pressure. Further, the material 
removal rate curves for different mixture ratios appear to 
intersect the nozzle pressure axis at around 55 MT/m 2 (gauge) . 

4.6.2 ITozzle Wear 

Hozzle wear profiles for split nozzles of specified 
geometrjr (Fig. 4*4) show that the wear in the convergent 
region of the nozzle is negligibly small compared to the 
straight region of the nozzle flow channel. The diameters 
of nozzle flow channel were measured along the length and 
diametral wear rates were evaluated. The diametral wear 
rate W^ can now be defined as, 

final diameter of nozzle _ previous diameter of 
flow channel (mm) “ nozzle flow channel (mm) 

d time (sec) 

(4.2) 


Variation of the diametral wear rate at different nozzle 
sections of straight region of the nozzle flow channel for 
different operating times are shown in Fig. 4.15. Figure 
shows that the diametral wear rate decreases with the increase 
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in operating time. It is further shown that the wear rate 
increases from the inlet towards the exit giving a maximum 
at the exit. 

The nozzle wear when expressed in terms of wear 
volume is thus a function of length and diameter of the 
nozzle. These parameters influence not only the wear volume 
lout also the wear rate since they affect the flow pattern 
of uhe particles inside the nozzle. In order to study the 
effect of the nozzle parameters on the wear rate, nozzle 
wear has "been plotted in terms of the dimensionless wear 
index W* which is defined as, 

final volume of nozzle ^ initial volume of nozzle 
_ ^low channel (nun ) _ flow channel (mm 3 ) 

initial volume of nozzle flow ch ann el 

(mm 3 ) 

(4.3) 

The wear index is plotted as a function of operating 
time for nozzles of various lengths in Fig. 4.16. The shape 
of the curve is similar for various nozzle lengths and indi- 
cates that the wear index initially increases linearly and 
then saturates at a specific value of operating time. As 
expected the wear index increases with increasing nozzle 
length (Fig. 4.17). 
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The variation of wear index with operating time 
for various nozzle diameters is shown in Pig. 4.18. The 
wear index for a specific operating time decreases appre- 
ciably when the nozzle diameter is increased (Pig. 4.19). 

Pig. 4.20 shows the variation of the wear index 
with operating time for various nozzle entrance angles. 

The initial wear rate (slope of the wear index curve) is 
again high hut tends to saturate rapidly. The wear rate 
increases with decreasing nozzle entrance angle. 

Pig. 4.21 shows the variation of the wear index 
with operating time for various mixture ratios. The wear 
index initially increases linearly with operating time hut 
saturates after a specific operating time. The wear index 
increases linearly with increasing mixture ratio hut the 
slope of the curve also increases with increase in operating 
time (Pig. 4.22). 

The effect of grain size on wear index is shown in 
Pig. 4*23- Tests with particles of grain sizes 30 and 38 
pm indicate that increased wear index is obtained with 
coarse grains. 
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r ig 4.1 Schematic layout of abrasive jet machining set-up. 




Fig. 4.2 (a) Experimental set-up of abrasive 
jet machining. 



Fig. 4.2(a) Photograph of mixing chamber, secondary abrasive 
chamber and vibrator unit . 
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Fig. 4.7 Machined cavity profile at different 

stand-off distances (a) 2mm, (b) 6mm, 
(c) 10mm, (d) I4mm,(e) 16 mm and 
(f ) 20 mm. d p : 30 pm ; CX : 0.148 ; nozzle 
pressure : 147.15 kN/m 2 (gauge) ; cutting 
time : 60 sec. 
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Fig. 4.16 Variation of the wear index with the operating 
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Variation of the wear index with the operating time 
for different nozzle exit diameters (a) 162 mm, 

(b) 1 17 mm ,(c) 0.9 mm and (d) 0.8mm Hardness 50 
HRc , In : 12 mm , : 30° d p : 38 pm , <*: 0.201. 



0.5 . 10 L5 2.0 

Nozzle diameter, mm 

Fig. 4.19 Variation of the wear index with the nozzle 
diameter for different operating times(a)5min 
(b) 10 min and (c) 15 min . Hardness : 50HRc * 
l n : 12 mm* £ : 30°$ dp : 38 pm* o( ; 0.201 . 
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CHAPTER 5 


DISCUSSION 


5.1 Erosion 

As mentioned earlier, material removal in AJM 
occurs due to erosive action of abrasive particles impact- 
ing at high, velocity onto the work surface. The mechanics 
of material removal in such situations is quite complex and 
difficult to analyse "because of many unknown and inter- 
dependent parameters. Prom one-dimensional steady state 
two-phase flow analysis presented in Chapter 2 the velocity 
of impacting particle can "be estimated. Considering quasi- 
static indentation under elastic-plastic response of target 
material and using semi-empirical fracture mechanics rela- 
tionships a theoretical model of erosion rate for "brittle 
materials such as glass, ceramics etc. has "been developed 
in Chapter 3. 

Equation (3.40) shows that "beside material parameters, 
the erosion rate in "brittle materials under solid particle 
impact is a strong function of the velocity of impacting 
particle. In most of the earlier investigations [25,46,50, 
63» 75 ], the emphasis was put on the velocity exponent 
controlling the erosion rate in "brittle materials. In the 
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present theory, the velocity exponent of 2.75 has been 
obtained. Earlier Weiderhom and Hockey [75] and Sheldon 
and Einnie [46] experimentally obtained the velocity 
exponent values as 2.5 and 5 , respectively for erosion of 
glass. Other brittle erosion theories [25, 65 ] give the 
velocity exponent in the range of 2.4 to 5.2. These 
theories, however, do not account for the effect of 
particle velocity threshold below which no erosion takes 
place. The present analysis includes this effect. 

Equation (5.44) gives the estimate of volumetric 

material removal rate in MM. It contains the term u 

pc 

which is essentially a material dependent parameter and 
can be evaluated from equation (5.8). To obtain the 
particle velocity at the jet centre line u at a specified 
stand-off distance, flow governing equations of gas- 
particle mixture flow in nozzle [equations (2.68), (2.69), 
(2.71), (2.75), (2.74), (2.75) and (2.84)] and in the jet 
region [equations (2.108), (2.110) and (2.112)] can be 
solved numerically taking particle velocity ratio c to be 
0.955 at the nozzle inlet as obtained in Appendix -VI. 
following the solution technique given in Chapter 2 the 
particle velocity at the jet centre line for various 
stand-off distances in a free air-abrasive jet has been 
computed for the nozzle geometry shown in Eig. 1.1 and 
for the input parameters listed in Appendix-VII. 
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Evaluation of volumetric material removal rate from 
equation (5.44) requires the knowledge of X . For the 
purpose of theoretical analysis* X was assumed to he a 
function of dimensionless stand-off distance x*^. "Values 
of X evaluated from the experimental data (Fig. 5.1) indi- 
cate that X can he expressed as 

X = xj (c 6 + c ? x* ) , (5.1) 

where 0 r and c_ are constants. Here dimensionless stand- 
6 ( 

off distance x* is given hy 



(5.2) 


where x, is the stand-off distance and d is the nozzle 
d no 

exit diameter. Constants Cg and c evaluated from the 

experimental data, Table 5*1» indicate that for the range 

of mixture ratio used the variations in their values are 

small. For the range of mixture ratio considered the 

-5 

average values of Cg and may he taken as 8.53x10 and 
-1.72xl0~ 4 . 


Using the proposed erosion model the variation of 
theoretical material removal rate with stand-off distance 
has been plotted in Fig. 5.2. For qualitative comparison 
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the experiment al data are also plotted which, show satis- 
factory agreement "between the theory and experiments. 
Experimental as well as theoretical results (Figs. 4.8 
and 5.2) show that the material removal rate first increa- 
ses with increase in stand-off distance and then decreases 
giving an optimum value. This pattern of curve is essen- 
tially due to the combination of inertia effects in sus- 
pended particles and particle-particle interaction effect . 
The carrier fluid attains a maximum velocity at the nozzle 
exit but abrasive particles being heavier tend to lag 
behind and attain the maximum velocity some distance away 
from the nozzle exit. On exit from the nozzle the jet 
flares up due to entrainment of the surrounding stagnant 
fluid and the carrier fluid velocity decreases. At some 
stage in the jet region the carrier fluid velocity and the 
velocity of abra,sive particles attain the same value i.e., 
the velocity of slip reduces to zero. Eeyond this the jet 
flaros significantly and the particle velocity decreases 
end eventually reaches a value (u . < u ) when no erosion 

JJLii 

can tahe place. At this stage the abrasive particles and 
worle material will essentially have elastic Impact. Equation 
(3.44) indicates that the erosion rate is a function of 
impingement velocity. Hence, the volume of material removed 
will increase as impingement velocity of abrasive particle 
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increases. With, the flaring of the jet, although the 
particle "velocity decreases, the effective number of 
impact per unit time increases due to decreased particle- 
particle interaction effect and the impingement takes 
place over a larger area. The combined effect of these 
parameters i.e. impingement velocity, particle-particle 
interaction effect and impingement area "will give an 
optimum value of stand-off distance for maximum material 
removal rate. 

Pigs. 4.8 and 5.2 further reveal that the maximum 
material removal rate at various values of mixture ratio 
appears to occur at around the same value of stand-off 
distance. Thus, the mixture ratio does not seem to have 
any effect on the value of stand-off distance correspond- 
ing to the maximum material removal rate. Althotigh the 
velocity of abrasive particle decreases with increase in 
mixture ratio at a specified stand-off distance, but the 
number of particles impacting the work surface increases, 
giving increased erosion rate. Further, in tho range of 
mixture ratios used it appears that the combined effect of 
particle-particle interaction and the velocity of abrasive 
particle do not change significantly since the maximum 
material removal rate at various values of mixture ratio 
are at almost the same stand-off distance. 
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The volume /weight of the eroded material is of 
importance in operations such as polishing, dehurring 
etc., while in micro-drilling and cutting operations it 
is the erosion depth (depth of penetration) which is 
more relevant. The variation of penetration rate with 
stand-off distance is shown in figs. 4-6 and 4.9. It 
is clear that the penetration rate also attains an optimum 
value hut the corresponding stand-off distance is much 
smaller than the value where maximum material removal 
rate is obtained. The penetration rate essentially 
depends on the impingement velocity and the particle- 
particle interaction effect, whereas material removal 
rate is also influenced by the area of impingement. The 
entrainment of the surrounding fluid results in flaring 
of the jet which causes further non-uniform velocity 
distribution of the solid particles across the jet sec- 
tion (Fig. 2.2). The particles at the centre will have 
higher velocity than those at the periphery of the jet. 

Due to this non-uniform velocity distribution of the 
solid particles in the jet the maximum value of the aver- 
age velocity of particles occurs at a smaller stand-off 
distance than that where the centre line velocity of the 
particle in the jet attains the maximum. The combined 
effect of particle-particle interaction and the maximum 
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value of the average velocity of the solid particles 
gives the optimum penetration rate. Fig. 4.6 indicates 
that the maximum penetration rate is attained at a stand™ 
off distance of 7.2 mm which corresponds to approximat ely 
10 times the nozzle exit diameter. Thus, this point -will 
correspond to the point where the combined effect of maxi- 
mum average velocity and effective number of impacts per 
unit time (inverse of particle-particle interaction effect) 
attain the maximum. Beyond this point the particle-particle 
interaction effect further decreases causing increased 
number of impacts per unit time but the decrease in average 
particle velocity is more, resulting in decreased penetra- 
tion rate. 

Fig. 4.9 shows that at low mixture ratio (0.095) 
the maximum penetration occurs at approximately 7.5 times 
the nozzle exit diameter, but when the mixture ratio is 
increased to 0.268 the maximum penetration occurs at 
approximately 14.5 times the nozzle exit diameter. This 
is because of the delayed acceleration of abrasive parti- 
cles when the mixture ratio is high. Flow analysis of 
gas-solid particle mixture jet (Sec. 2.7.2 and 2.7.3) 
reveals that the distance at which the particle in the 
jet attains the maximum velocity depends on the length 
of the potential core (x ) of the jet. Equation (2.97) 
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shows that the potential core length in gas-solid particle 
mixture jet depends on the mixture ratio, higher mixture 
ratio gives larger potential core length. Thus, the dis- 
tance at which the average value of particle velocity in 
the jet attains the maximum will depend on the potential 
core length and thereby on the mixture ratio. The pene- 
tration rate depends on the average value of the particle 
velocity at any section of the air-abrasive jet and there- 
fore at higher values of mixture ratio maximum penetration 
rate will occur at larger stand-off distance. The increase 
in penetration rate with increasing mixture ratio is dtie 
to the increase in number of particles taking part in ero- 
sion in a specified time. 

The top diameter of eroded cavity increases with 
increase in mixture ratio as well as stand-off distance 
as shown in Fig. 4.10. This effect is duo to the flaring 
of the jet which increases with increase in stand-off 
distance as well as with increase in mixture ratio. 

Fig. 4.11 indicates that coarse particles give 
higher values of volume erosion factor than fine particles. 
Equation (5.44) indicates that for a specified mass flow 
rate of abrasive particles the volumetric material removal 
rate is proportional to the impacting particle diameter 
and thus coarse particles will give higher volume erosion 
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factor. For a specified mixture ratio the number of 
particles per unit volume of mixture will he less for 
the coarse particle resulting in decreased particle- 
particle interaction compared to fine particle. Coarse 
particle although attains low velocity hut due to dec- 
reased particle-particle interaction effect it gives 
more number of impacts onto the work surface in a speci- 
fied time giving higher material removal rate. The 
velocity achieved hy fine particles is higher hut increa- 
sed particle-particle interaction effect results in lower 
volume erosion factor. It appears from Fig. 4.11 that in 
the size range of particles used, the combined effect of 
particle size and number of impacts onto the work surface 
in a specified time results in increased volume erosion 
factor with coarse grains. The net effect of particle- 
particle interaction, velocity and the diameter of particle 
also appears to he such that the optimum stand-off distance 
for maximum value of volume erosion factor remains more or 
less unaffected. 

Fig. 4.12 shows that the optimum material removal 
rate increases almost linearly with increase in mixture 
ratio. The higher the mixture ratio, the greater will he 

the number of particles per unit volume of the mixture. 

rj 

Thus, the number of particles impacting the work surface 
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increases with increase in mixture ratio. The proportio- 
nality relationship in Pig. 4.12 will, however, exist as 
long as the particles remain in suspension with the carrier 
fluid and the particle-particle interaction effect is small. 
Thus, for low mixture ratio values the material removal rate 
can he taken to he directly proportional to the mixture 
ratio. At higher values of mixture ratio, the curve will 
have a drooping characteristic. This is because at higher 
mixture ratios the quantity of the carrier fluid will not 
he sufficient to impart the required momentum and energy 
to the abrasive particle and the particle may tend to 
settle. Also, at higher mixture ratios the particle- 
particle interaction effect increases and the particle 
velocity decreases since less energy is available from 
the carrier fluid. These factors will cause decreased 
material removal rate at higher mixture ratios. Por the 
same mixture ratio coarse grains give higher material 
removal rates than fine grains because of larger impacting 
particles and decreased particle-particle interaction 
effect. Also the proportionality limit of mixture ratio 
for fino particles is more than that for the coarse parti- 
cles due to the inertia effects. Por any specified mixture 
ratio, the coarse particles being heavier have larger 
velocity slip resulting in lower particle velocity. Thus, 
the limiting value of mixture ratio upto which the material 
removal rate is proportional is lower for coarse particle. 
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Fig. 4.13 shows the variation of penetration rate 
with stand-off distance for different particle sizes. The 
coarse particles being heavier accelerate slower and hence, 
the penetration rate attains a maximum at a larger value 
of stand-off distance. Further, the coarse particles have 
small particle-particle interaction effect thus have more 
number of impacts on the workpiece surface in a given time. 
The depth of erosion per particle impact is proportional 
to the lateral crack depth produced under the impact which 
in turn is proportional to the impacting particle diameter 
[equation (3.29)]. Thus, coarse particles will give higher 
penetration rate. 

The effect of nozzle pressure on volumetric material 
removal rate under specified machining conditions is shown 
in Fig. 4.14. It can be seen that as pressure increases 
the rate of material removal increases but it tends to 
saturate beyond a certain value. In convergent -straight 
nozzles the maximum velocity attainable at nozzle exit is 
sonic velocity [119]. Thus, increasing the pressure above 
a critical value would not increase the jet velocity beyond 
the sonic velocity and the material removal rate will tend 
to saturate. The critical pressure limit , however, depends 
on the mixture ratio, high mixture ratio requiring higher 
pressure to reach the sonic velocity [90]. Increased mixture 
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increases the number of particles exposed to acceleration 
by the carrier fluid resulting in increased particle drag. 
Also increased number of particles imposes a proportionally 
higher frictional drop on air flow through nozzle. These 
effects will cause the pressure required for achieving a 
given velocity at the nozzle exit to increase with increas- 
ing mixture ratio as seen in fig. 4.14. In order to achieve 
erosion, a minimum impact velocity (u ) is required. This 

threshold velocity (u ) will vary with the properties of 

pc 

the impacting particle and the workpiece material [equation 
(3.8)1. fig. 4.14 indicates that the threshold velocity is 
achieved when nozzle pressure is around 55 Ml/m 2 (gauge). 

5.2 Uozzle Wear 

Ductile erosion studies [12-18] have revealed that 
the erosion rate is a function of the velocity and the 
angle of impingement of the impacting particles [equation 
(1.3) ]. Experiments using split nozzles have revealed that 
most of the wear occurs in the straight passage at the 
nozzle exit. The wear is negligibly small in the convergent 
region of the nozzle, further, the wear rate increases from 
the entry to the exit of the straight region of the nozzle 
as shown in fig. 4.4. Theoretical estimate of solid particle 
velocity inside the nozzle of convergent -straight geometry as 
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shown in Pig. 4.4(a), is obtained using the procedure 
given in Chapter 2 for the input conditions listed in 
Appendix-VII and is plotted in Pig. 5.3 as a function 
of nozzle length to exit diameter ratio. Pig. 5.3 shows 
that the increase in the velocity of solid particle is 
negligibly small in the convergent region of the nozzle 
whereas faster increase in velocity is seen in the 
straight region of the nozzle. Purther, the solid parti- 
cle velocity increases with increase in length to exit 
diameter ratio of the nozzle. This is because of the 
inertia effect of the susp ended particles. Since the 
carrier fluid is lighter it accelerates faster and attains 
its maximum velocity at the end of convergent region of 
the nozzle, whereas heavier solid particles tend to lag 
behind the fluid and continue to accelerate even in the 
straight region of the nozzle. In the initial part of 
the straight region the solid particle accelerates faster 
because of large velocity of slip, whereas in later part 
the acceleration rate decreases and attains a steady 
state value. As the wear volume depends strongly on the 
velocity, it increases as the particle moves in the straight 
region of the nozzle. Pigs. 4.16 and 4.17 provide further 
evidence that the wear rate is higher for longer nozzle 
lengths . 
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The particle velocity within the nozzle is also a 
function of the nozzle diameter. The variation of theore- 
tical velocity of particle with the length of nozzle for 
various nozzle diameters is shown in Pig. 5.4. For a 
specified inlet diameter and entrance angle of the conver- 
gent part s the length of the convergent region of the 
nozzle will he greater for smaller exit diameters. For 
particulate flow the solid particles accelerate because 
of the exchange of momentum and energy due to interaction 
with the carrier fluid. Because of the effects of inertia 
a finite length is required to accelerate the solid parti- 
cles. Thus, the particle velocity will he greater for 
smaller exit diameters and vice-versa. The nozzle diameter 
increases during use, resulting in a decreased particle 
velocity which in turn causes a decrease in the wear rate. 
This is clearly shown in Figs. 4.16 and 4.18 where the 
curves show a decrease in the wear rate (slope of the wear 
index curve) with time. The wear index is initially almost 
linear and then saturates. This is because the velocity 
of the particle decreases as nozzle diameter increases dtie 
to erosion and hence the wear rate decreasc-s from its ini- 
tial high value. The effect of nozzle diameter on wear 
index is shown in Fig. 4.19- Since the particle velocity 
is greater for smaller exit diameter (Fig. 5-4), the result- 
ing wear rate is high. 
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Pig. 4.20 shows the effect of nozzle entrance angle 
on the wear index. The wear increases with decreasing 
entrance angle because the effective length of the conver- 
gent region of the nozzle increases resulting in increased 
particle velocity. Thus the wear rate must increase with 
decreasing entrance angle. 

The mixture ratio is another important parameter 
affecting erosion in gas--solid flow. Pig. 4.22 shows that 
the wear rate increases when the mixture ratio is increased. 
This is obvious from the fact that the number of particle 
impacting on the nozzle surface will increase with increas- 
ing mixture ratio and hence the erosion rate will increase. 
The erosion rate, however, is unlikely to increase conti- 
nuously since the carrier fluid at higher mixture ratios, 
may not be able to impart sufficient momentum to accelerate 
the solid particles. Thus, acceleration of solid particles 
beyond a specific mixture ratio will be low, resulting in 
a low particle velocity and a decreased wear rate. 

The effect of grain size on nozzle wear has also been 
investigated by comparing the results obtained using abra- 
sive particles of 30 and 33 pm. The erosion of the material 
is a function not only of the impingement velocity but also 
of the size of the impinging particles because the volume of 
the cavity eroded by a single particle is a function of the 
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contact diameter between the impacting particle and the 
target surface [24]. The contact diameter, in turn, 
is a function of the particle diameter, the properties of 
the particle and the target materials and the impingement 
velocity [108]. A larger particle will remove more mate- 
rial than a smaller particle at the same impact velocity. 
The inertia of a large particle is however greater and 
hence the velocity attained will be lower. The number of 
particles impinging on the surface at a given mixture ratio 
also depends on the particle size. Thus, smaller particles 
will undergo more impacts but , since the contact diameter 
is small, they will erode a smaller volume, larger grains, 
in contrast, will -undergo fewer impacts but, because of 
their larger contact diameter, will erode a larger volume. 
The effect of particle size on erosion rate is therefore 
a combined function of particle velocity, particle size 
and the number of impacting particles. In the range of 
particle size used, larger particles appear to give higher 


erosion rate 
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TABLE 5.1 


Values of the Experimental Constants Cg and c ^ 


Mixture 

ratio 

0.095 

0.148 

0.201 

0.268 


c 6 


Crj 


8.26 x 10~ 3 
8.63 x 10~ 3 
9.16 x 10“ 3 
8.07 x 10“ 3 


- 1.81 x 10 ^ 

- 1.80 x 10" 4 

- 1.93 x 10~ 4 
~ 1.34 x 10~ 4 


Average values ; 

c g = 8.53 x 10" 3 
= -1.72 x 10“ 4 



X*x ri {05 





10 15 20 25 30 


)imensionfess stand-off distance , 


ariation of X with dimensionless stand-off 
stance for different mixture ratios(a)0.095» 

) 0.148, (c) 0.201, (d) 0.268 and(e)mean curve 
* BO um : nozzle oressure ; 147.15 kN/m 2 
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CHAPTER 6 


CONCLUSIONS AND SCOPE POR FURTHER 
RESEARCH 


6.1 Conclusions 

liras ive jet machining can be used for polishing, 
de burring and other finishing operations where the rate 
of material removal is important and for operations such 
as micro- -drilling and cutting where the erosion depth is 
more relevant. The theoretical analysis for brittle ero- 
sion predicts that erosion rate depends on density, size 
and velocity of impacting particles as well as hardness, 
fracture toughness and shear modulus of the workpiece 
material. The velocity of impacting particles, however, 
is the main parameter governing the erosion rate. The 
value of velocity exponent in the erosion equation is 2.75 
for typical glass and ceramic materials. There also exists 
a threshold velocity below which no erosion occurs and its 
value depends on workpiece hardness, density of impacting 
particle and elastic constants of workpiece and particle 
materials. Further, the theoretical analysis of two-phase 
flow through nozzle as well as jet reveals that the velo- 
city of solid particle in air-abrasive jet depends on 
particle size, density of particle material, mixture ratio 
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and carrier gas pressure. Solid particles being heavier 
in inertia continue to accelerate even after the nozzle 
exit and attain the maximum velocity at some distance away 
from the nozzle exit. The analysis therefore indicates 
that the erosion rate depends on stand-off distance, mixture 
ratio, particle size and carrier gas pressure. 

The experimental study of effects of these parameters 
on volumetric material removal rate as well as penetration 
rate have been successfully carried out on the experimental 
abrasive jet machine specially designed and fabricated for 
the purpose. Results indicate that both material removal 
rate and penetration rate first increase with increase in 
stand-off distance and then decrease giving optimum values. 
Theoretical curves also indicate this and there appears 
qualitative agreement between the two. Maximum values of 
material removal rate and penetration rate are, however, 
obtained at different values of stand-off distance. Increase 
in mixture ratio causes increase in material removal rate 
as well as penetration rate. In the range of mixture ratios 
used during the experiments, the material removal rate 
appears to vary linearly with mixture ratio. This linearity 
limit appears to vary with grain size. The maximum value 
of material removal rate appears to occur at the same value 
of stand-off distance, while the maximum penetration rate 



186 


occurs at different values of stand-off distance when 
mixture ratio is varied. The size of abrasive particle 
affects the value of stand-off distance for maximum 
penetration rate, but does not seem to affect the value 
of stand-off distance for maximum material removal rate. 

It appears from the result that operations such as deburr- 
ing and polishing should be carried out at larger stand- 
off distances, whereas smaller stand-off distances should 
be txsed for cutting, micro-drilling etc. 

Experiments performed using split convergent- 
straight nozzles clearly indicate that nozzle wear increa- 
ses almost linearly from entry to exit. The profile of 
the straight flow passage after wear closely approximates 
a truncated cone. The wear volume could therefore be 
estimated with reasonable accuracy from measurements of 
the entrance and exit diameters of the straight passage 
of the nozzle. Wear in the convergent portion appears to 
be negligibly small. Experiments clearly indicate that 
nozzle wear increases with increasing nozzle length and 
decreasing nozzle diameter. Similarly, nozzle wear increases 
when the entrance angle of the nozzle is decreased. Nozzle 
wear rate also increases with increase in particle size and 


mixture ratio. 
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6.2 Scope for Further Research 

The mechanics of material removal in AJM and the 
wear characteristics of the nozzle used in the process 
are not fully understood and the present work can he exten- 
ded for the case of axially symmetric three dimensional 
analysis of two-phase flow. The axi-symmetric impinging 
jet analysis can he carried out to account for particle- 
particle interaction and velocity, and effective number of 
impacts can he obtained at different stand-off distances 
for realistic estimate of material removal rate as well 
as penetration rate. 

The work can he ex-bended to study the tapering 
effect due to flaring of the jet, an inherent character- 
istic of the process. This effect can perhaps he minimized 
by using co-axial twin jet of air and air-abrasive mixture. 
The invest igat ion can he carried out theoretically as well 
as experimentally to minimize jet flaring. 

Practically nothing appears to he known about 
nozzle design. The nozzle used imparts high velocity to 
the abrasive particles and thus the nozzle geometry has 
significant influence on the pex’formance characteristic of 
the process. The work can he extended to investigate the 
effects of nozzle parameters on process performance and 


nozzle life 
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APPENDIX - I 

D IMSF S 10 1'TLE S S FORM OF THE CROSS-- SECTIONAL AREA OF 

NOZZLE AND ITS DERIVATIVE 


Considering Fig. 1.1 the distances and 
measured along the nozzle length from inlet to the end 
of I convergent region and to the end of II convergent 
region can be expressed as 


and 


(A 


nl 


d n2^ 


2 tan jB^ 


1. = In 


(a. 


n2 


a ) 

no 


2 tan p, 


( 1 . 1 ) 


( 1 . 2 ) 


respectively. Here d ^ 1 , d 2 and d nQ are the diameters 
at the inlet, at the end of I convergent region and at 
the nozzle exit, respectively. |3^ and arG scm i~ 
entrance angle for the I and II convergent regions of the 

- V 

nozzle . 

For any section at a distance x from the nozzle inlet 
in I convergent region, the diameter of nozzle passage d 
can be express od. as 

d = d - 2 x tan p n for x £ 1-, . (1.3) 

n nl 1 i 
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The area A of the nozzle at the considered section is 
n 

given as 


A . = ( d n -j - 2 x tanjB^) 2 for x < 1^ (1.4) 

Using dimensionless variables discussed in Chapter 2 
[Sec. 2.5.2 equation (2.64)], equations (I.l), (1.2) 
and (1.4) can be written in dimensionless form as 


and 


* _ ^ d nl " d n2^ 

L 1 ' 2 tanp^ 


( 1 . 5 ) 


.* = d al , (tm h ~ * an ^ d n2 

l 2 • 2tanp 1 2 tanp^ tan|3 2 


2tan^ 2 

( 1 . 6 ) 


^ = (i- 


2 x tanp 1 




) 2 for z* < 1* . (1.7) 


Here dimensionless diameters g. and 'g ^ can be expressed 
as 


a 


-)r 


nl 


nl 


no 


( 1 * 8 ) 
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and 




(1.9) 


Differentiation of equation (1.7) with, respect to f 
yields 




= -4(1 - 


2x* tan(3- 


nl 


0 (- 


tanP- 


* 

l nl 


■) for x <l n 


(I. 10) 


Similarly in the II convergent region the diameter d^ of 
the nozzle passage can be expressed as 


tanpp 

4 n = <d nl " + d n2 " 

for 1 1 < x <l 2 (1. 11) 


The area A^ can be expressed as 


A 


71 


n 


[(a. 


t an(3 

+ & n2 ~ 2xtan £ol 


hi n2 t anp- 


for 1 ± < x < 1 2 (1.12) 


Using dimensionless variables of equation (2.64) equation 
(1.12) can be written as 
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* 

K 


a 


= [(1 


* 

n2 

* 

nl 


tanp, 

tanp- 


* 


+ (■ 


l n2, 2x tan ^2 l2 

f a* 

nl nl 


for 1* < x*<l* (1.13) 

Differentiating equation (1.13) with respect to x* gives 



a 


* 

X 


4 tanj3 0 d 0 tanp 0 

2 r / x _n2\ _2 

u - 1 tanp. 


Co 2x*tanp 0 

+ (-&)- — ~2] 


nl 


nl 


nl 


nl 


for 1* < x* < 1* (1.14) 

After II convergent region i.e., in the straight 
region of the flow passage the nozzle diameter is constant 
and is equal to the exit diameter of the nozzle d nQ . In 
this region the nozzle flow area in dimensionless form can 
he expressed as 




for x* > 1* 


(1.15) 


Differentiation of equation (1.15) with respect to f 
yields 


* 



for x* > l| 


o 


(1.16) 
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APPENDIX - II 

MAXIMUM CONTACT PONCE DURING- ELASTIC COLLISION 

Using Hertz relation [108] the contact force P 

during elastic collision of a spherical particle with a 

flat body can be related to the elastic deformation depth 

z , impacting particle diameter d and elastic constants 
e p 

of the impacting particle and the target material as 


E 


e 




■)] (' 


2 ) 4 - 


1- V- 


P 



(ii.i) 


During collision when the impacting particle comes 
in contact with the flat body, the compressive force P g is 
set up which resists the penetration of the impacting 
particle into the target surface. The equation of motion 
of the colliding particle can be written as 


and the term 


( ? d3 
v 6 p 


) _ 
j dt 

- P 

e 

(II. 2) 

du p 

dt"' as 





li 

Pd 

to 

d_ 

dt 

dz„ 

.(-—£) - 
l dt ; 

at 2 

(II. 3) 


i 
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Equations (II- l), (II. 2) and (II. 3) give 


d 2 Z 


at 


2 r 4 / 
" z e [ 3 ( : 


1- 


a 


)](^)V(f dp P"p) . 


"P. 


E 


(II. 4) 


dz @ 

Mult iplication of on both sides of equation (II. 4) 

gives 


dz 

* ^ 


4 '? 

71 '■ 


1- 

(t 


v. 


1- v‘ 


+ _*) (p d^) 


5 

a 2 


2 

-] z 2 dz 


P 


P P 


(II. 5) 


When impacting particle touches the target surface, the 
deformation rate will be equal to the velocity of impinge- 
ment and at full penetration the velocity becomes zero. 
Thus the above equation (II. 5) can be integrated as 


S 

u 


dz 

* 4(av> 2 


P 


4 f 2 

TC 



1 


^ Z e^max 


i — v* 2 
+ t^)(p d 2 ) 

v 1 x p p ' 


] / 


z 


3/2 



(II. 6) 
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Here (z )„ rt _ is the maximum value of elastic deformation depth 
at the end of impact cycle when particle velocity reaches 
a zero value. Integration of equation (II. 6) gives 


^ z c ^max 


5n P T 


L l6fZ 




(II. 7) 


Substitution of z o = (z ) in eq_uation (II. l) gives 

c? msi X 


(F ) 
e'max 


,20% - w - 

( 'T ■°p ) S'- 
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1" t 

1 92( _t.,. + 


1- at 


2 6 

t5 5 
■'J < • 


P- 


P 


(II. 8) 
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APPENDIX - III 


SAIT AT ION VELOCITY POP HORIZONTAL CONVEYING- 


The saltation velocity u for horizontal conveying 

& s 

of solid particles can be estimated using the procedure 

presented by Zenz [115]. Prom Pig. III.l, the minimum 

velocity (u ) . needed to convey a single particle of 
gs mm 

largest size is obtained. Por this the constants A^ and 
A^ are evaluated using equations 


lf.° ikJVfaal 3 5, cm/s eo 


3 p 


gs 


(III.l) 


and 


L 4 


= [ 


3 n; 


4g p ( p - p ) 
°c gs v p gs 


P . 


cm 


(III. 2) 


Here \x a is the viscosity of the carrier fluid g/cm-sec, 

‘o 

g is the acceleration due to gravity cm/sec 2 , P ^ is the 
material density of the solid particles g/cm 3 and p is 

o o 

the material density of the carrier fluid g/cm 3 . 

In the present case aluminium oxide particles (max. 
50 pm) are to be transported using air at the temperature 
of 30° and pressure of 0.5 MN/m 2 (abs) as carrier fluid. 
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For this condition the values cf different variables used 
in equation (III.l) and (III. 2) are 

jly = 1.86 z 10~ 4 g/cm~sec, 

P gs = 57.44 z 10" 4 g/cm 3 , 

P = 3.85 g/cm 3 , 

Jr 

and g = 981 cm/sec 2 . 

Using equations (III.l) and (III. 2) the constants and 

A , for the above conditions are 
4 

A„ = 30.797 cm/sec 

5 

and 

A^ = 1.0516 z 10-3 cm. 

For d = 0.005 cm , 



£.005 

o To 010 5 16 


4.755 


and from Fig. III.l for spherical particles 



min 



4.76 
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which gives (u ) . = 146.57 cm/sec. Also from Pig. III.l 

the slope of the tangent (n^ )on the considered curve at this 

point is zero. Since the data of Pig. III.l pertains to 

the flow in 6. 35 cm tube, the correct saltation velocity 

(u )'. for the actual size of the tube can be obtained 
v gs mm 

from the following equation : 


^ U gs 'rain 


^ u gs ^min 


( 6.35 ) 


0.4 


(III. 3) 


where d, is the actual size of the tube, 
t 

In the pi’esent case the minimum diameter of the 
nozzle is 0.7 mm, the system inlet pressure and temperature 
are 0.5 MM/m 2 (abs) and 30°C , the nozzle exit pressure and 
temperature are 0.101 MF/m 2 (abs) and 30°C , respectively. 

The carrier fluid velocity at the nozzle exit in the present 
case is taken to be 34000 cm/sec (see Sec. 4.3)* Ike conti- 
nuity relation for the flow between nozzle exit and the tube 
gives 


d 


t 


34000 (d^ Q p g0 ) ^ 

p ( u ) 

gs gs min 


(III. 4) 


where d is the nozzle exit diameter in cm, p is the 
"no S° 

density of carrier fluid in g/cm 3 at nozzle exit conditions 

and p is the density of carrier fluid in g/cm 3 at system 

gs 

working conditions in the tube. 
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Substit uf ion of numerical values of o , d . p „ 

• go no 'gs 

and (u ) . into equation (III. 4) gives 
gs mm B 

P 34000 x (0.07) 2 x I.l64xl0" 3 
d , = [ — ~ — ~ — — - — j 2 = 0.48 cm 

x 5. 744x10“ 3 x 146.57 

Prom equation (III* 3) the corrected value of minimum 
saltation velocity for the tube diameter (d^) of 0.48 cm 
is 


^ u gs ''min 


146.57 ( 


0.4S nQ « 4 

6 . 35 


52.17 cm/sec • 


This value of saltation velocity (u._) pertains 

gp m ixi 

to single particle conveying which is further corrected for 
the actual mixture ratio through equation 


M P , p gs 


(f) u 

P 

P 


gs 


A tess_ 

5 (u ) . 

v gs'mm 


( u gs ^mm/ 

I ~ 


(III. 5) 


where is the mass flow rate of solid particle in g/sec 
and M is the mass flow rate of carrier fluid in g/sec in 

O 

the mixture. The constant is 

A.- = 21.4 n**'* for n- > 0.068 

5 3 j 


= 0.32 


for -0.11 < n^ 1 -0.068 


(III. 6) 
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For the value of mixture ratio considered (a = 0.4) » 

the ratio M comes out to be 0.667. Since n, = o, the 
P g 3 

constant may be ta^en to be 0.32. Substitution of these 
values in equation (III, 5) gives 


0 . 667 ( 51 - 4 ^^) ^ . 

3.965 gs 


(u - 52.17) 

°. 32 


which after solution gives u = 61.93 cm/sec. Thus, in 

gs 

the present case to transport 50 pm diameter aluminium 
oxide spherical particles of uniform size through a hori- 
zontal tube at mixture ratio ( cc ) of 0.4, the saltation 
velocity comes out to be 0.62 m/sec. 
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APPENDIX - IV 

DESEN OP MIXING- CHAMBER AND MIXTURE CONVEYING- TUBE 

The present set-up has been designed for the system 

inlet pressure of 0.5 MN/m 2 (abs) and temperature of 3 >0°C 

-.3 

at which the carrier fluid (air) density p = 5.744xXO 

o ^ 

g/cm 3 . The exit pressure and temperature are 0.101 MhT/m 2 
(abs) and 30°C i respect iveljr at which carrier fluid density 

p = 1. l64xl0” J g/cm 3 . The minimum size of the nozzl e for 

3° 

which the set-up is to be designed is 0.07 cm. As men-fioned 
in Sec. 4.3 the exit velocity of the carrier fluid for* the 
above stated conditions of flow is 34000 cm /sec at the nozzle 
exit . The minimum safe velocity of the carrier gas in. the 
mixing chamber (u^.^) and the mixture conveying tube 
are 62 cm/sec and 362 cm/sec, respectively (see Sec. 4-3). 

The continuity relation for the flow between nozzle 
exit and mixing chamber gives 

a _ [ 3 4000 dnoj jgo ^ ^ (IT .13 

“ " ^3 u geh 

where d is the diameter of the mixing chamber in cm. 
m _ 

Substitution of the numerical values of d no » Pgo* ^gs 
u v, in equation (IV.l) gives 
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d _ r 34000 x (0.07) 2 x 1 . 164x10'" 3 
m " 5 • 744x10” 3 x 62 

= 0.7379 cm. 

The safe value for the mixing chamber diameter is 7.3 mm. 

Similarly the continuity relation for the flow 
between mixture conveying tube and nozzle exit gives 


a 


mt 


= c 


34000 


1 2 P~ 
no go 




U 


gs gst 


(IV. 2) 


where d^ is the mixture conveying tube diameter in cm. 

Subs it ut ion of numerical values of the variables in eq.ua- 
t ion ( IV . 2 ) gives 


d _ j- 34000 x (0 . 07) 2 x 1.164x10“° 
mt 5. 744x10" 3 x 362 

= 0.3054 cm. 

Thus, the safe value for the mixture conveying tube diameter 
is 3 mm. 
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APPENDIX - V 


CALCULATION OF MIXTURE RATIO 


The manufacturer's calibration of the flow meter 
for air at 30°C temperature indicates 1 division (1 mm) 
to be 50.39 ml /min at system pressure of 372.454 kU/m 2 
(gauge) and 51.92 ml /min at system pressure of 344.863 
k?T/m 2 (gauge). The air density p is given by 


o 


0 . 001293 
(1 + 0.00567 T) 



(Y.l) 


whore P is the pressure in terms of cm of mercury and 
T is the temperature in °C. For 30°C temperature equation 
(7.1) gives p . = 51.265 x 10“ 4 g/ml at a pressure of 
34 4.863 kU/m 2 (gauge) and = 54.435 x 10" 4 g/ml at a 

o 

pressure of 372.453 kU/m 2 (gauge). 

Prom manufacturer's calibration chart supplied with 
the meter the equivalent number of divisions (mm) are 
obtained for the flow meter reading (reading at the centre 
of flow meter ball) which for appropriate values of system 
pressure with corresponding values of density gives the air 
flow rate in g/min. Using equation (l.l) the mixture ratio 
can be obtained for the known mass flow rate of abrasive 
particles . 
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APPENDIX - VI 

PARTICLE VELOCITY RATIO AT NOZZLE INLET 

Usirg equation (2.2) and (2.3) the motion of 
spherical particle (following the particle path) in a 
uniform stream of gas can he expressed as 


dUj 

dt 


P 


5 ( P -f ) ( r- > <" g - V 2 

Pp s p 


4 


(VI. 1) 


The solution of equation (VI. l) for particle 

initially at rest (x=o, t =o, u =o) is 

P P 




s 


(VI. 2) 


g P 


where 


t _ 


2 ( A ) °d 

4 v " ' d 

p p p 


(VI. 3 ) 


and x is the distance required to accelerate the particle 
from zero velocity to a velocity u . 

ir 

Equation (VI. 2) can be rewritten as 




Uy 


U 


u = 


^ • 


(VI. 4) 



216 


* 

In the present study for = 4.973 x 10 kg /sec 
a = 0.148, T f/1 = 303°E, P ]L = 24.8536 x 10 4 E/m 2 (abs). 


-A 


[a. = 0.186 x 10 ‘ E b / m‘ 

& 


’ d nl 


-5 

3 x 10 m, equation 


(2.66) gives u ^ = 14.17 m/sec and equation (2.67) gives 
P = 2.858 kg/m 3 . For aluminium oxide particles with 

£> -L 

material density p = 3850 kg/n 3 aid 30 pm size, equation (2.5 

P 

gives the particle Reynold 1 s number = 65.32. Equations 
(2.6) and (2.7) give the corrected drag coefficient 0^=1. 28. 
Substitution of this value of G^ in equation (VI. 3), gives 


= 23.75. 


The tube length for conveying air-abrasive mixture 

from mixing chamber to nozzle inlet is 0.5 m (i.e. x=0.5 m). 

Thus, for x = 0.5 m and § = 23-75, equation (VI. 4) gives 

the particle velocity ratio c as 

s 


c 


s 





0.935 . 


(VI, 5) 



APPENDIX ~ VII 


INPUT PARAMETERS POR THEORETICAL EVALUATION 
OE PARTICLE VELOCITY AMD EROSION RATE 

Carrier Eluid 

Carrier gas s Air 

Gas constant (R ) : 287 J/kg°E 

O 

Nozzle pressure (P^) ; 248. 475 kU /m 2 (abs) 

Inlet temperature (T ): 303°K 

gJ- 

Viscosity (p cj .) : 1.86 x 10“^ M s/m 2 

g 

Thermal conductivity (kg): 26.753 x 10*"*^ J /m 

Specific heat (c ) : 1004.6 J /kg°K 

i TO 

Density at nozzle exit (p ) : 1.164 kg/m 3 

o ^ 

Prandtl number (Sp r ) ! 0.701 

Abrasive Particles 

Type : Aluminium oxide 

Modiolus of elasticity (E ) : 36.55 x 10"^ 3f/m 2 
Poisson’s ratio (o> ) : 0.25 

Material density (P ) : 3850 kg/m 3 

Specific heat (c ) : 795.5 J/kg K 

Particle temperature (T^): 303°E 
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Work Material 


Type : 

Modulus of elasticity (E^ ) : 
Fracture toughness (K ) : 

L/ 

Vicker’s hardness (H^ ) : 

Shear modxilus (G , ) : 

Poisson's ratio ( ) ; 


Glass 

6.9 z 10 10 E/m 2 

7.5 x 10 5 F/m 5 / 2 

5.6 x 10 9 F/m 2 
2.76 x 10 10 F/m 2 
0.25 


S.Fo. Mass flow rate of particles Mixture ratio 

( kg/sec ) 


1 . 

3.193 


io -5 

0.095 

2 . 

4.974 

X 

IQ" 5 

0.148 

3. 

6.755 

X 

10 " 5 

0.201 

4. 

9.000 

X 

10 " 5 

0.268 
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